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This  thesis  introduces  a  climate  index  based  on  radiative  transfer 
theory  and  derived  from  the  spectral  radiances  typically  used  to  re¬ 
trieve  tenperature  profiles.  It  is  assumed  that  clouds  and  climate  are 
closely  related  and  a  change  in  one  will  result  in  a  change  in  the 
other.  Because  the  index  developed  in  this  paper  is  primarily  a  func¬ 
tion  of  the  cloud  distribution  and  also  dependent  upon  temperature  and 
moisture  distributions,  it  may  be  used  as  a  climate  index.  The  advan¬ 
tage  is  that  the  index  is  more  accurately  retrieved  from  satellite 
data  than  cloudiness  per  se.  This  index,  hereafter  referred  to  as  the 
VIRES  index  ( for  Vertical  Infrared  Radiative  Emitting  Structure ),  is 
based  upon  the  shape  and  relative  magnitude  of  the  broadband  weighting 
function  of  the  infrared  radiative  transfer  equation.  The  broadband 
weighting  curves  (VIRES  curves)  are  retrieved  from  simulated  satellite 
infrared  sounder  data  (spectral  radiances).  The  VIRES  index  is  a  use¬ 
ful  means  of  classifying  the  VIRES  curves.  This  thesis  describes  the 
retrieval  procedure  and  investigates  error  sensitivities  of  this  method. 
Advantages  of  the  approach  include;  day  and  night  capability,  no  need 
to  know  cloud  radiative  properties,  retrieval  ability  when  cloud  frac¬ 
tion  or  cloud  emittance  is  less  than  1.0,  retrieved  information  below 
sensor  resolution,  minimum  influence  on  the  index  from  low  tropospheric 
retrieval  errors,  and  the  ability  to  statistically  discriminate  between 
VIRES  curves  retrieved  from  different  scenes.  Operational  VIRES  index 
retrieval  strategies  and  a  number  of  specific  applications  are  also 
proposed. 
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ABSTRACT  OF  THESIS 


A  CLIMATE  INDEX  DERIVED  FROM  SATELLITE  MEASURED 
SPECTRAL  INFRARED  RADIATION 

This  thesis  introduces  a  climate  index  based  on  radiative  transfer 
theory  and  derived  from  the  spectral  radiances  typically  used  to  re¬ 
trieve  tengerature  profiles.  It  is  assumed  that  clouds  and  climate  are 
closely  related  and  a  change  in  one  will  result  in  a  change  in  the 

other.  Because  the  index  developed  in  this  paper  is  a  function  of  the 

cloud,  temperature,  and  moisture  distributions,  it  may  be  used  as  a 
climate  index.  The  advantage  is  that  the  index  is  more  accurately  re¬ 
trieved  from  satellite  data  than  cloudiness  per  se.  This  index,  here¬ 

after  referred  to  as  the  VIRES  index  (for  Vertical  Infrared  Radiative 
Emitting  Structure),  is  based  upon  the  shape  and  relative  magnitude  of 
the  broadband  weighting  function  of  the  infrared  radiative  transfer 
equation.  The  broadband  weighting  curves  are  retrieved  from  simulated 
satellite  Infrared  sounder  data  (spectral  radiances).  This  thesis 
describes  the  retrieval  procedure  and  investigates  error  sensitivities 
of  this  method.  It  also  proposes  index  measuring  options  and  possible 
applications  of  the  VIRES  index. 

Results  indicate  that  the  VIRES  approach  is  a  very  effective  use 
of  satellite  radiometer  measurements.  Retrieval  advantages  include; 
day  and  night  capability,  no  need  to  know  cloud  radiative  properties, 
retrieval  ability  when  cloud  fraction  or  cloud  emittance  is  less  than 
1.0,  minimal  geometric  assumptions,  retrieved  information  below  sensor 
resolution,  and  minimum  influence  on  the  index  from  low  tropospheric 
retrieval  errors.  These  advantages  along  with  the  approach  of  compos¬ 
iting  scenes  for  an  average  VIRES  curve  greatly  reduce  the  retrieval 


ill 


sensitivity  to  the  kind  of  errors  found  under  assumed  normal  operating 
conditions.  A  detailed  error  analysis  indicated  that  the  most  impor¬ 
tant  error  sources  are  instrument  system  noise,  and  ill  specified 
temperature  and  humidity  profiles.  Accurate  VIRES  retrievals  are  il¬ 
lustrated  under  a  number  of  different  error  and  atmospheric  conditions. 
Furthermore,  a  statistical  technique  used  to  successfully  discriminate 
between  VIRES  curves  derived  for  different  atmospheric  conditions  is 
discussed. 

Operational  VIRES  index  retrieval  strategies  and  a  number  of 
specific  applications  are  proposed.  It  is  suggested  that  the  index  be 
derived  from  geostationary  satellite  data  and  averaged  to  provide  week¬ 
ly  regional  values.  These  index  values  would  be  used  in  a  regional 
climate  monitoring  mode.  They  would  also  be  useful  for  verification  of 
climate  model  generated  infrared  radiation  to  space  values. 

Michael  David  Abel,  Captain  USAF 
Atmospheric  Science  Department 
Colorado  State  University 
Fort  Collins,  Colorado  80523 
Spring,  1981 
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I.  INTRODUCTION 

This  paper  introduces  and  describes  a  climate  index  called  the 
VIRES  index.  VIRES  (pronounced  vi*res)  is  the  acronym  for  Vertical 
Infrared  Radiative  Emitting  Structure  and  quite  by  chance  it  is  also 
the  plural  form  of  the  Latin  word  vis,  meaning  forces  or  powers.  Thus, 
the  acronym  seems  especially  appropriate  since  the  atmosphere's  VIRES 
is  one  of  the  major  forcing  factors  behind  the  earth’s  climate.  The 
VIRES  index  is  related  to  the  earth's  climate  through  the  earth's  radi¬ 
ation  budget  and  therefore,  can  be  considered  a  climate  index.  The 
logic  of  this  may  be  stated  as  follows.  The  climate  system  is  deter¬ 
mined  by  the  energy  input  to  the  system  and  the  distribution,  trans¬ 
formation,  and  storage  of  energy  in  various  forms  within  the  system. 
These  processes  are  mirrored  in  the  components  of  the  earth's  radiation 
budget,  one  of  which  is  the  outgoing  emitted  thermal  radiation  (COSPAR 
Report  to  ICSU  and  JOC,  1978b).  This  cooling  to  space  is  described  by 
the  VIRES  which  is  primarily  a  function  of  cloud  distribution. 

Atmospheric  observation  has  been  and  continues  to  be  central  to 
the  progress  of  atmospheric  science.  Better  observations  remain  one  of 
the  needs  of  the  discipline.  The  recent  introduction  of  meteorological 
satellite  systems  has  contributed  significantly  to  the  growing  need  of 
monitoring  world-wide  weather  variables.  Satellites  not  only  have 
world-wide  coverage  capability  with  good  horizontal  and  time  resolu¬ 
tion,  they  also  have  a  second  advantage.  Large  numbers  of  observations 
are  made  with  the  same  instrument  increasing  the  integrity  and  compar¬ 
ability  of  such  measurements  (Houghton,  1979).  The  relative  if  not  the 
absolute  accuracy  of  the  satellite  measurement  is  high.  More  attention 
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to  absolute  accuracy  has  been  paid  to  the  sounding  type  instruments 
than  to  any  other  satellite  instrument)  (i.e.  ITPR,  SIRS,  IRIS). 

The  scientific  community  is  still  seeking  the  best  ways  to  use  and 
apply  satellite  data.  Most  of  the  recent  successes  in  using  satellite 
data  in  a  quantitative  global  way  in  the  atmospheric  sciences  have  come 
in  the  areas  of  solar  constant  measurement  and  radiation  budget  measure¬ 
ments  (Heath,  1973;  Smith  et  al.  1977;  Vonder  Haar  and  Oort,  1973).  Ap¬ 
plication  of  satellite  data  to  these  problems  is  fairly  straightforward 
since  the  principal  satellite  instrument  is  a  radiometer  and  the  meas¬ 
urement  is  a  spectral  or  broadband  irradiance.  Other  important  areas 
of  research  include  inference  of  temperature  and  humidity  profiles,  as 
well  as  cloud  and  wind  determinations  from  satellite  radiance  values. 

The  suitability  of  satellite  observations  is  reduced  since  the  required 
information  must  be  inferred  from  the  radiance  values  measured  remotely 
at  the  satellite  and  from  the  appropriate  geometric  and  radiative  trans¬ 
fer  principles. 

The  objective  of  this  study  is  to  describe  the  Vertical  Infrared 
Radiative  Emitting  Structure  (VIRES)  of  the  atmosphere  by  using  simu¬ 
lated  satellite  spectral  radiation  measurements.  A  process  is  describ¬ 
ed  that  accomplishes  this  objective  using  specified  cloud  radiative 
properties  and  mean  temperature  and  gaseous  atmospheric  profiles. 

These  findings  are  used  to  examine  the  feasibility  of  using  infrared 
radiative  transfer  weighting  curves  (which  describe  the  atmosphere's 
VIRES  by  defining  how  the  atmosphere  cools  to  space)  as  a  climate  in¬ 
dex.  This  index  would  be  principally  dependent  on  climatological 
cloudiness,  and  its  variability  could  be  regarded  as  an  indicator  of 
climate  variance.  The  strength  of  this  approach  is  that  it  deals 


directly  with  the  radiative  aspect  of  the  problem  thus  circumventing 
the  need  to  infer  specific  individual  clouds  directly  using  the  typical 
parameters  of  height,  base,  and  amount.  The  distinction  of  this  tech¬ 
nique  is  that  computationally  one  may  be  able  to  simply,  accurately, 
and  quickly  archive  the  atmospheric  VIRES  as  a  manifestation  of  global 
cloudiness  in  a  form  that  is  useful  for  monitoring  climate  change,  or 
for  validating  the  statistical  characteristics  of  cooling  to  space 
confuted  by  climate  models.  Many  factors  account  for  the  strong  po¬ 
tential  of  this  approach. 

Day  to  day  variations  as  well  as  longer  period  variations  in  the 
atmospheres  VIRES  (which  is  reflected  in  satellite  measured  earth 
radiances)  are  primarily  caused  by  clouds.  Many  techniques  using 
satellite  data,  some  of  which  are  discussed  in  the  next  section,  have 
been  developed  to  infer  cloudiness  in  the  standard  sense.  Of  course 
any  inference  of  clouds  using  such  data  will  by  definition  be  a  kind  of 
radiative  measure  of  cloudiness  with  the  drawback  that  specific  radia¬ 
tive  cloud  properties  must  be  assumed  before  results  are  possible.  Of 
course  it  is  desirable  to  make  as  few  a  priori  assumptions  as  possible 
when  analyzing  the  data  for  the  purpose  of  obtaining  reliable  cloud  in¬ 
formation.  By  using  a  radiance  measure  of  cloudiness  as  proposed  in 
this  paper  one  increases  the  compatibility  between  the  satellite  meas¬ 
urement  and  the  quantity  labeled  cloudiness.  Furthermore,  by  using  a 
unique  form  of  the  technique  commonly  called  the  infrared  sounder  cloud 
retrieval  method  for  a  single  field  of  view,  we  minimize  the  assump¬ 
tions  about  the  spatial  scales  and  geometry,  and  about  the  cloud  radia¬ 
tive  properties  while  avoiding  many  time  consuming  calculations  in¬ 
volving  iterations  through  the  radiative  transfer  equation.  This  new 
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method  is  used  to  find  two  curve  shape  parameters  which  define  the 
atmosphere's  VIRES  in  terms  of  the  Infrared  weighting  curve. 

Following  the  chapter  on  background  information,  the  specifics  of 
the  approach  are  described.  For  orientation  purposes  Figures  la,b  are 
provided.  As  described  in  Figure  la  the  procedure  starts  with  satel¬ 
lite  spectral  data  (in  this  case  simulated  data)  in  the  15  Urn  CO^  ab¬ 
sorption  bands  and  the  10-12  ym  window  band.  These  radiances  contain 
Information  on  the  VIRES  of  the  atmosphere.  Computationally,  relative 
importance  is  placed  on  the  radiances  depending  upon  where  in  the 
vertical  most  of  its  energy  originates.  By  assuming  known  or  measured 
gaseous  and  temperature  profiles  one  may  interpret  the  scene  radiative- 
ly  by  solving  for  the  weighting  function  peak  due  to  radiatively  spec¬ 
ified  effective  clouds  and  the  fractional  weighting  of  a  totally  over¬ 
cast  effective  cloud  scene  versus  a  totally  clear  scene.  This  scheme 
results  in  a  spectral  weighting  curve  shape  specified  by  the  two  vari¬ 
ables  mentioned  in  Figure  la.  By  design  the  curve  shape  is  not  depen¬ 
dent  on  the  cloud  radiative  properties  specified.  For  example,  if  a 
cloud  covering  the  entire  satellite-sensed  scene  is  specified  as  opaque 
(black)  when  its  emittance  is  only  0.5,  the  routine  will  compute  a 
proper  weighting  function  peak  height  (Pwf)  due  to  the  cloud  with  0.5 
fractional  weighting  (a) .  These  two  curve  shape  variables  will  give 
the  correct  weighting  function  curve,  the  same  curve  one  gets  from  an 
overcast  case  and  cloud  emittance  of  0.5.  Of  course,  if  one  insists 
on  interpreting  the  weighting  function  peak  and  fractional  weighting 
as  cloud  top  height  and  cloud  fraction,  the  accuracy  of  the  cloud  frac¬ 
tion  value  is  strongly  dependent  on  how  close  the  assigned  cloud  emit¬ 
tance  is  to  the  true  cloud  emittance. 
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Figure  la.  Schematic  diagram  showing  computational  steps  leading  to  the  proposed  atmospheric 
VIBES  climate  index. 


Pressure  (mb)  Pressure  (mb) 


Figure  lb.  Schematic  diagram  showing  the  relationship  between  four 

atmospheric  cloud  scenarios  and  the  resulting  VIRES  curves 
and  VIRES  indices  (reference  Chapter  VI-B). 
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Once  the  two  variables  determining  weighting  curve  shape  are 
found  they  may  be  used  in  a  broadband  transfer  equation  to  give  a 
broadband  weighting  curve  (see  Figure  lb).  This  curve  combined  with 
the  respective  temperature  profile  describes  the  vertical  structure  of 
atmospheric  cooling  to  space.  As  a  test,  the  value  of  the  earth's 
emittance  to  space  calculated  from  the  derived  weighting  curve  may  be 
compared  to  a  satellite  measured  value  of  the  same  quantity.  As  part 
of  this  study  an  error  analysis  is  done  to  evaluate  the  influence  of 
several  assumptions  on  the  results.  In  addition,  the  important  ques¬ 
tion  of  time  and  space  averaging  of  the  quantities  discussed  above  is 
also  addressed.  A  technique  for  discriminating  between  weighting 
function  curves,  and  a  related  climate  index  is  discussed.  Finally, 
the  strengths  and  limitations  of  using  IR  broadband  weighting  curves 
as  a  climatic  index  and  representation  of  cloudiness  will  be  examined. 


II .  BACKGROUND  INFORMATION 


There  have  been  many  attempts  to  deduce  cloud  cover,  structure, 
and  radiative  properties  from  satellite  data.  Table  1  is  an  outline  of 
most  of  the  approaches,  all  of  which  seek  to  define  inferred  cloudiness 
in  standard  terms.  A  short  summary  of  the  more  notable  research  fol¬ 
lows.  However,  a  critique  of  specific  approach  shortcomings  is  not 
attempted. 

Using  visible  wavelength  values  Miller  and  Feddes  (1971),  have 
related  brightness  measurements  to  cloud  amount.  Analysis  of  cloud 
amount  and  type  from  satellite  pictures  (nephanalysis)  has  been  done 
using  the  'eyeball*  method  (Clapp,  1964).  Another  more  objective 
method  combines  pictures  from  two  geostationary  satellites  to  give  a 
steroscopic  view  and  measure  of  cloud  height  (Dalton  et  al.  1979).  Of 
course  these  methods  are  limited  to  daylight  observations  and  are  de¬ 
graded  by  variable  and  cloud  look-a-like  surface  reflectance. 

Infra-red  window  data  combined  with  simplifying  assumptions 
(which  greatly  reduce  the  accuracy  and  applicability)  and  a  'known' 
temperature  profile  when  used  with  appropriate  radiation  laws  will  give 
estimates  of  cloud  height  or  cloud  fraction  (Koffler  et  al.  1973).  One 
also  needs  to  know  or  estimate  cloud  radiative  characteristics.  An 
example  of  a  technique  that  uses  this  type  of  data  is  the  adjacent 
field  of  view  method  described  by  Smith  et  al.  (1970).  They  use  the 
derived  cloud  information  to  construct  clear  column  radiance  profiles 
as  part  of  a  temperature  profile  retrieval  process.  Exact  cloud  loca¬ 
tion  is  still  temperature  profile  dependent.  Another  totally  different 
approach  is  described  by  Rao  (1970).  He  statistically  relates 


METHODS  OF  OBSERVING  CLOUDS  FROM  SATELLITES 

I.  Visible  Wavelength  Observations 

A.  Reflected  Solar  Radiances 

1.  Threshold 

2.  Weighted  histogram 

B .  Nephanalys is 

C.  Stereoscopic 

II.  Infra-red  Wavelength  Observations 

A.  Window  Measurements 

1.  Threshold 

2.  TD 

B.  CO^  Gas  Band  Measurements 

1.  Single  field  of  view  -  RTE  iterations 

a.  Radiance  ratioing 

b .  Minimization 

2.  Single  field  of  view  -  empirical  RTE* 

III.  Vis/IR  Combined  Observations 

A.  Dual  channel 

B.  Bi-spectral 

C.  2-D  histogram 

♦Developed  and  employed  in  this  paper 


Table  1.  An  outline  of  satellite  cloud  retrieval  methods. 
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radiative  values  to  differences  between  surface  and  cloud  top  tempera¬ 
tures  (TD)  over  ocean  areas.  However,  the  most  straightforward  ap¬ 
proach  is  to  simply  relate  cloud  top  temperature  and  a  known  or  assumed 
temperature  profile  (assuming  one  measures  a  single  overcast  cloud 
layer) .  An  example  of  this  threshold  approach  is  the  work  of  Cox  and 
Griffith  (1978)  using  GATE  data.  More  recently  Canpbell  et  al.  (1980) 
have  used  geosynchronous  satellite  IR  window  observations  to  produce 
area  cloud  top  distribution  profiles.  They  wish  to  assess  the  impact 
of  the  diurnal  and  spatial  changes  of  these  distributions  on  the  earth- 
atmosphere  radiation  budget. 

Other  methods  seek  to  improve  accuracy  by  combining  solar  bright¬ 
ness  and  IR  window  information.  However,  while  improving  the  accuracy 
over  taking  each  method  separately,  one  must  be  content  with  the  limi¬ 
tations  of  both  methods.  A  good  example  of  this  is  the  bi-spectral 
technique  of  Reynolds  and  Vonder  Haar  (1977)  and  expanded  upon  by 
Mendola  and  Cox  (1978) .  They  solve  simultaneously  a  set  of  budget  type 
radiative  equations .  They  also  use  a  method  described  by  Shenk  and 
Curran  (1973)  to  improve  retrieval  of  cirrus  clouds.  Other  methods  are 
often  referred  to  as  dual  channel.  For  example,  in  another  paper 
Reynolds  et  al.  (1978)  describe  a  technique  for  discriminating  differ¬ 
ent  cloud  types  by  visible  and  IR  image  subtraction.  A  similar  ap¬ 
proach,  called  the  2-D  histogram  method,  is  outlined  by  Smith  (1978). 
This  last  technique  is  designed  to  handle  large  quantities  of  data  very 
quickly.  Of  course  result  accuracy  is  sacrificed  for  speed.  However, 
for  climatological  applications  the  results  may  be  useful. 

Because  this  paper  presents  a  technique  using  COj  gas  band  meas¬ 
urements,  this  approach  will  be  discussed  in  greater  detail  below. 
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First  it  should  be  mentioned  that  there  are  other  techniques  for  deter¬ 
mining  cloud  characteristics  from  satellite  measurements  that  do  not 
fit  neatly  into  the  categories  outlined  in  Table  1.  For  example,  it 
may  be  possible  in  the  future  to  use  Ildar  techniques  from  space  like 
those  described  by  Platt  (1979)  for  ground  based  units.  Microwave 
measurements  may  be  used  with  SW  and/or  IR  measurements  in  a  tri-channel 
or  dual  channel  approach  (Yeh  and  Llou,  1980) .  Following  this  reasoning 
a  3-D  histogram  technique  has  been  suggested.  A  method  using  IR  window 
and  water  vapor  channels  is  being  pursued  by  Chen  et  al.  (1980).  Fi¬ 
nally,  a  technique  that  uses  spectral  infrared  measurements  from  limb 
scanning  is  described  by  Taylor  (1974)  and  by  Remsberg  et  al.  (1980). 

To  conclude  this  section  the  single  field  of  view  CO2  gas  band 
measurement  technique  will  now  be  discussed.  It  is  also  referred  to  as 
the  infrared  sounder  cloud  retrieval  method.  This  approach  has  many 
advantages.  It  requires  the  fewest  a  priori  assumptions  while  provid¬ 
ing  day  and  night  capability.  However,  it  does  have  problems  detecting 
low  clouds.  More  details  on  the  assumptions  involved  and  limitations 
will  be  given  later. 

The  IR  RTE  in  integro-differential  form  represents  the  backbone 
of  this  method  and  is  given  below. 

lnp  (1) 

o 

L(V,0)  =  e8(v)  B(V,  Ts)  T (v,  e,  p8)  +  f  B(v,  T[p])  aT-^p  9-,-P)  d  lnp 

lnP8 

2  -1 

where  L  is  spectral  radiance  in  W/m  sr  cm  ,6  is  solar  zenith  angle 

and  e  is  surface  spectral  emittance  usually  taken  to  be  1  for  10-15 
s 

ym  wavelengths .  The  Planck  Function  B  is  given  below. 


B(V,  T) 


V  /[exp(c2v/T)-l] 


where  V  is  wavenumber  in  cm  \  T  is  temperature  in  aK,  and  and  c2 
are  constants.  The  equation  for  transmittance  T  is  given  below. 


t(v,  0,  p) 


K  (v,p)  sec  6  dp] 


where  q  is  the  gas  mass  mixing  ratio,  g  is  the  acceleration  due  to 
gravity,  p  is  pressure  with  pq  being  pressure  at  the  top  of  the  atmos¬ 
phere  and  p  being  pressure  at  the  surface,  and  K  is  the  gas  absorption 
s 

3t 

coefficient.  Also,  g-"^p  is  commonly  referred  to  as  the  weighting 
function.  Equation  1  may  be  rewritten  as  Eq.  (2)  for  the  case  of 
opaque  clouds  with  cloud  top  height  at  p£  for  cloud  fraction  N,  and 
1-N  clear  sky  (cs) . 


lnpo 

L(V,  6)  -  N  Jb(v,  Tc)  T(v,  0,  pc)  +  f  B(v,  T[p])  |-^  d  lnpj 

lnpc  ' 


+ 


(1-N) 


|b(v,  Ts)  t(v,  0, 


inpo 

ps)  +  J  B(V, 
lnPs 


3  T 
3  lnp 


d  lnp) 


LVcld  + 


(1-N)  L. 


VC3 


(2) 


Two  basic  techniques  using  these  equations  to  solve  for  cloud  proper¬ 
ties  have  been  proposed.  One  is  known  as  the  radiance  ratioing  method. 
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It  is  described  in  Smith  and  Woolf  (1976),  McCleese  and  Wilson  (1976), 
and  Smith  and  Platt  (1978).  It  is  also  used  by  Wlelicki  and  Coakley 
(1980),  who  have  described  its  applicability  and  limitations  in  detail. 
An  outline  of  this  method  follows. 

First  rewrite  Eq.  (2)  for  grey  clouds  using  the  relationship  o  * 
e  N  for  a  spectral  radiance  of  wavenumber  Vl. 


\l  -  aLvldd  +  «-“>  lvlc.  (3) 

where  a  and  Lylcld  are  the  unknowns  and  Lvlcld  depends  only  on  pc.  Re¬ 
arrange  Eq.  (3). 


(Lvl  ’  Lvlcs)  '  “  (Lvlcld  “  KjIcs^*  (4 

To  have  one  equation  with  one  unknown  (Pc) ,  ratio  Eq.  (4)  for  two  dif¬ 
ferent  wavenunber  radiances  and  eliminate  a  assuming  N  ■  N 


(L. 


'vl 


(L 


V2 


LylC8 

Lv2cs 


) 

) 


(LVlcld  ~  Lylcs  ^ 
<Lv2cld  -  Ntfcs* 


(5) 


Iterate  through  different  pc’s  until  the  and  that  satisfy  Eq. 
(5)  are  found.  Finally  take  the  just  found  and  solve  Eq.  (4)  for 


a. 
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A  second  method  Is  described  by  Chahine  (1974).  It  Is  based  upon 
minimization  of  the  RMS  difference  between  the  observed  radiances  and 
calculated  radiances  that  are  a  function  of  cloud  top  pressure  and  ef¬ 
fective  cloud  fraction.  An  iterative  scheme  is  used  to  choose  the 
cloud  pressure  and  fraction  used  in  the  radiative  transfer  equation 
calculations.  While  the  first  technique  is  designed  to  use  only  two 
band  radiances,  the  second  method  may  use  more  than  two  bands.  Requi¬ 
site  lengthy  radiative  transfer  calculations  are  a  disadvantage  to 
operational  use  of  this  method. 

Both  of  these  techniques  and  the  one  developed  in  this  paper  as¬ 
sume  the  profile  T(p)  is  known  which  Implies  the  clear  sky  spectral 

radiance  L  is  known.  All  three  methods  assume  spectral  band  emit- 
Vcs 

tances  (e^)  are  equal,  thus  the  relationship  is  assumed  true. 

All  three  assume  the  clear  sky  spectral  band  weighting  functions  are 
known  and  are  not  identical  to  each  other.  All  three  methods  assume 
the  satellite  radiance  measured  comes  from  a  scene  that  contains  only 
the  representative  grey  body  cloud  top  pressure.  Measurements  in 
either  the  4.3  ym  or  15  ym  CO^  band  can  be  used.  Sometimes  the  window 
channel  (11  ym)  is  also  used  with  the  CO2  bands  even  though  they  are 
widely  separated  in  wavenumber  (McCleese  and  Wilson,  1976).  In  this 
case,  only  as  e^-KL  does  evl  -  (Yamamoto  et  al.  1970).  For  this 
reason  poor  results  can  be  expected  using  the  window  and  COj  channels 
together  to  detect  nonblack  clouds.  Notice  that  all  terminology  used 
in  this  paper  is  consistent  with  recommendations  of  the  IAMAP  Radiation 
Commission  except  that  V  is  used  for  wavenumber  instead  of  K  (Raschke, 
1978).  The  next  chapter  describes  the  third  technique  mentioned  above. 
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III.  SHAPE  PARAMETER  RETRIEVAL  THEORY 

As  stated  earlier,  the  basic  approach  used  in  this  paper  to  deter¬ 
mine  the  weighting  function  curve  shape  parameters  is  commonly  called 
the  infrared  sounder  cloud  retrieval  method  for  a  single  field  of  view. 
The  commonly  used  procedures  are  described  in  detail  by  Chahlne  (1975) 
and  by  Smith  and  Platt  (1978).  Chapter  II  of  this  paper  contains  a 
brief  review  of  their  techniques.  Because  of  the  limited  number  of 
assumptions  needed  and  the  day-night  capability,  the  CO2  band  approach 
is  without  question  the  most  accurate  for  determining  high  and  middle 
cloud  information.  Low  cloud  Information  in  tropical  atmospheres  is 
limited  primarily  due  to  the  high  concentration  and  emittance  of  water 
vapor  in  the  lower  troposphere  which  masks  the  clouds'  radiant  signa¬ 
ture.  However,  as  will  be  shown  this  limitation  is  minimized  by  re¬ 
trieving  the  atmospheres  VIRES  instead  of  clouds  per  se.  A  brief 
description  of  the  radiative  theory  behind  this  method  follows. 

Equation  1  gives  the  IR  integral  form  of  the  RTE  for  a  plane- 
parallel  clear  sky  atmosphere  with  no  scatter  under  the  assumption  of 
local  thermodynamic  equilibrium.  A  plot  of  the  variation  of  the  trans¬ 
mittance  with  respect  to  pressure  is  called  the  weighting  function. 

For  radiation  measured  in  the  absorption  band  or  window  band  in  the 
absence  of  clouds,  the  transmittance  Is  a  known  function  (with  slight 
dependence  on  temperature  and  water  vapor  profiles)  as  is  the  weight¬ 
ing  function. 

Equation  2  represents  the  satellite-received  spectral  radiance 
[L(v,9)]  from  a  scene  containing  N  fraction  of  radiatively  black  clouds 
For  nonblack  clouds  N  takes  on  a  different  meaning  and  is  replaced  by 


a,  an  effective  cloud  fraction  (see  Eq.  3).  Because  reflectance  by 
black  or  nonblack  clouds  is  assumed  to  be  small  compared  to  emittance 
in  the  IR  bands  of  concern,  effective  emittance  e(v)  and  transmittance 
T (v)  are  related  by  e(v)  +  x(v)  ■  1  for  these  wavenumber  bands.  Con¬ 
sequently,  we  may  as  stated  above  replace  N  with  an  equivalent  frac¬ 
tional  cloud  cover  value  equal  to  the  product  e(v)  N  ■  a.  Thus,  true 
fractional  cloud  cover  cannot  be  derived  unless  the  cloud  emittance  is 
known. 

Neglecting  term  1^ ,  which  is  usually  calculated  from  assumed 

temperature  and  gas  profiles,  Eq.  (2)  has  unknowns;  p£  and  N.  If  we 

assume  e(v^)  *  e(v2)...  for  the  wavebands  of  interest,  then  a  =  e(v)  N 

will  be  the  same  value  in  Eq.  (3)  written  for  each  of  the  spectral 

radiances.  We  therefore,  have  a  system  of  at  least  two  equations  with 

two  unknowns,  pc  and  a.  So  far  we  have  assumed  knowledge  of  gaseous 

and  temperature  profiles  which  allow  us  to  compute  the  appropriate 

clear  sky  spectral  radiance  and  the  many  (p^)  values  (using 

specified  cloud  radiative  characteristics)  used  to  find  p  and  a.  We 

c 

also  assume  that  the  clouds  in  the  scene  all  have  approximately  the 
same  cloud  top  pressure  level  pc-  However,  we  do  not  need  to  assume 
the  scene  is  overcast,  nor  is  it  necessary  to  make  assumptions  about 
adjacent  scenes. 

At  this  point  the  method  used  in  this  paper  diverges  from  the  so- 

called  ratio  method  and  minimization  method  described  in  Chapter  II. 

Instead  of  using  Eq.  (3)  that  gives  satellite  measured  radiance  in 

terms  of  a,  L  ,  and  L,  . an  empirical  equation  with  L  in  terms  of 
Vcs  veld  V 

a,  L  ,  and  p  is  derived.  This  step  eliminates  the  need  to  solve  the 
vcs  *c 

IR  -  RTE  for  the  iteration  value  of  L,  ...  The  rationale  is  outlined 
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below  using  simulated  satellite  radiances.  The  spectral  radiative 
transfer  routine  and  specific  wavenumber  bands  and  cloud  radiative 
characteristics  used  are  described  in  Chapter  IV. 

The  emphasis  is  on  the  relationship  between  specified  effective 
clouds  and  the  weighting  function.  With  specified  effective  clouds 
present  the  weighting  function  is  drastically  changed.  For  example. 
Figure  2  shows  the  effect  graphically  using  wavenumber  747.5  cm  *  (5 
cm  *  band  width)  and  different  levels  of  overcast  cloudiness  for  a 
typical  tropical  atmosphere.  Cloud  emlttance  is  near  1  or  is  unity 
since  thick  clouds  are  specified  using  an  emlttance  model  (described 
later)  related  to  specified  cloud  water  content.  Figure  3  shows  the 
effect  of  effective  clouds  on  the  weighting  function  for  different 
amounts  of  cloudiness  for  a  tropical  atmosphere.  This  is  the  same 
effect  as  changing  cloud  emlttance  to  less  than  unity  in  an  overcast 
case.  In  other  words  a  is  the  important  shape  parameter  (a  =  Ne^) . 
Broadband  weighting  functions  show  the  same  general  characteristic  ef¬ 
fects  of  clouds,  except  for  low  altitude  effective  clouds.  Near  the 
ground,  water  vapor  (see  Figure  4)  acts  much  like  a  low  effective  cloud 
as  far  as  the  atmosphere's  VIRES  is  concerned.  This  is  particularly 
true  in  the  tropics.  Figures  3  and  5  show  that  small  a's  rdsult  in 
small  changes  in  the  weighting  curve  shape. 

For  each  wavenumber  interval  in  the  CC^  band  there  is  a  different 
shaped  weighting  function.  Wavenumbers  close  to  the  center  of  the  band 
show  clear  sky  weighting  function  peaks  near  the  tropopause  due  to 
stronger  line  absorption.  Wavenumbers  further  from  center  such  as 
747.5  cm-1  show  peaks  at  lower  levels  in  the  atmosphere.  When  there  is 
an  effective  cloud  present,  primary  or  secondary  peaks  (Pw^)  occur  near 
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the  effective  cloud  top  level  as  seen  in  Figures  2  and  3.  The  exact 
location  of  the  peak  is  a  function  of  model  vertical  resolution  and  in 
this  20  mb  vertical  resolution  model  it  is  found  10  mb  below  specified 
cloud  top  (pc  *  pwf  -  10). 

The  most  important  radiative  property  in  IR  bands  is  cloud  emit- 
tance  as  a  function  of  depth  into  the  cloud.  Investigation  shows  that 
for  a  model  vertical  resolution  of  20  n£*,  the  level  of  the  weighting 
function  peak  due  to  the  cloud  is  not  sensitive  to  cloud  emittance 
specifications.  Even  for  transluscent  clouds  the  weighting  function 
peak  (Pwf)  is  10  ni>  below  specified  effective  cloud  top  (pc) . 

To  find  an  empirical  relationship  between  L  s>  Pw£»  and  ct,  values 
of  weighting  function  peak  heights  (Pwf)  due  to  specified  effective 
clouds  were  plotted  against  the  corresponding  values  of  L  (pc>  a  =  1) 
for  a  given  temperature  humidity  profile.  Figure  6  shows  these  plots 
for  a  mid-latitude  summer  profile.  Figure  7a, b  contains  similar  plots 
for  a  tropical  atmosphere.  The  following  relationship  was  found: 

<t>  <nr>  ’  ttv  -  K*>  <6> 


The  constants  C  and  D  depend  on  spectral  wavenumber  and  atmospheric 
profile  used  although  for  small  profile  changes  (1  to  2°C  or  20-30% 
water  vapor)  they  are  nearly  constant.  The  method  used  to  derive  and 
apply  Eq.  (6)  is  the  topic  of  Chapter  IV. 


Empirical  RTE  curves  for  the  mid-latitude  summer  model  for  a  *  1  case. 
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IV.  SHAPE  PARAMETER  RETRIEVAL  PROCEDURE 

Radiances  received  at  satellite  level  are  simulated  using  a  spec¬ 
tral  radiative  transfer  equation  (RTE)  for  infrared  radiation  developed 
by  Cox  et  al.  (1976).  Spectral  absorption  data  are  taken  from  Elasser 
and  Culbertson  (1960),  Smith  (1969)  and  Bignell  (1970).  For  broadband 
infrared  calculations  a  routine  is  used  that  is  a  broadband  approxima¬ 
tion  to  a  rigorous  line  by  line  spectral  radiative  transfer  equation, 
and  which  is  described  by  Cox  et  al.  (1976),  and  by  Griffith  and  Cox 
(1977).  Both  sets  of  computer  code  were  modified  to  produce  the  out¬ 
put  requirements  of  this  research.  By  design  both  radiative  models 
are  computationally  fast  with  the  consequence  that  approximations  re¬ 
sult  in  decreased  accuracy.  For  example,  N20  and  CH^  absorption  is 
ignored  (Gupta  et  al.  1978).  Nevertheless,  the  principals  of  the 
method  described  below  are  not  dependent  on  the  absolute  accuracy  of 
the  radiative  calculations. 

Spectral  bands  chosen  for  use  in  this  research  are  typical  of 
those  used  on  the  satellite-borne  radiometers  called  VTPR  -  Vertical 
Temperature  Profile  Radiometer  (NOAA  2-5) ,  HIRS  -  High  Resolution  In¬ 
frared  Sounder  (Nimbus  6,  TIROS-N)  and  VAS-VISSR  Atmospheric  Sounder. 
Detailed  descriptions  of  the  instruments  may  be  found  in  McMillin 
et  al.  (1973),  Sissala  (1975),  and  Schwalb  (1978).  Table  2  gives  the 
central  wavenumber  of  the  bands  used  in  this  study.  These  represent 
typical  values  and  are  not  necessarily  the  optimum  ones.  A  5  cm  ^ 
band  width  is  used.  Table  2  also  gives  clear  sky  atmosphere  weighting 
function  properties  of  these  channels. 
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The  tropical  and  mid-latitude  atmospheric  variables  used  in  this 
study  are  given  in  Tables  3  and  4  respectively.  These  values  were 
taken  from  the  Handbook  of  Geophysics  and  Space  Environments  (1965)  and 
from  U.S.  Standard  Atmosphere  1976  (1976).  As  Tables  3  and  4  indicate, 
the  radiative  transfer  routine  is  run  with  a  20  mb  tropospheric  resolu¬ 
tion. 

Besides  gas  and  temperature  profiles,  cloud  radiative  properties 
must  be  specified.  As  stated  earlier  and  as  examined  in  more  detail 
later,  these  specified  cloud  characteristics  are  not  critical  to  ob¬ 
taining  the  proper  weighting  function  curve  shape  parameters  (p^f  and 
a) .  Figure  8  describes  the  emittance  model  used  in  this  research. 

Table  5  gives  the  specified  cloud  parameters.  Notice  that  each  of  the 
45  modeled  effective  clouds  is  100  mb  thick  (where  possible)  and  that 
below  300  mb  the  emittance  is  unity  (black  radiating  surfaces).  The 
effective  cloud  tops  range  from  100  mb  to  980  mb  at  20  mb  intervals. 

Since  all  parameters  have  been  defined,  radiance  values  can  now 
be  calculated  for  each  of  the  six  bands  given  in  Table  2.  Satellite 
received  radiance  values  are  simulated  for  each  atmosphere  and  wave 
band  for  forty-five  overcast  cases  and  one  clear  sky  case.  From  these 
230  spectral  radiance  values,  any  simulated  set  of  satellite  values 
for  a  given  atmospheric  profile  can  be  generated  using  Eq.  (3).  Fig¬ 
ures  6  and  7a  give  plots  of  Lv  vs.  for  the  overcast  case  (a  values 
given  in  Table  5)  simulated  in  this  way.  The  procedure  for  obtaining 
the  weighting  curve  shape  parameters  will  be  discussed  next. 

A  standard  least  square  linear  regression  model  (Snedecor  and 
Cochran,  1967)  is  applied  to  the  data  in  Figures  6  and  7a  using  a  log 
pressure  transformation.  The  basic  model  is  y  ■  b  exp  (mx)  and  in  its 
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pressure 

mb 

TEMPERATURE 

°K 

W  H20 
G/KG 

W  03 
OG/G 

W  C02 
G/KG 

1.6 

265.0 

.010 

6.200 

.486 

20.0 

225.6 

.020 

10.342 

.486 

40.0 

215.3 

.010 

4.429 

.486 

60.0 

205.6 

.010 

1.723 

.486 

80.0 

198.7 

.010 

.623 

.486 

100.0 

195.7 

.010 

.348 

.486 

120.0 

200.0 

.010 

.227 

.486 

140.0 

206.0 

.010 

.200 

.486 

160.0 

211.0 

.010 

.175 

.486 

180.0 

216.4 

.010 

.152 

.486 

200.0 

221.0 

.010 

.138 

.486 

220.0 

225.2 

.020 

.126 

.486 

240.0 

228.7 

.040 

.114 

.486 

260.0 

232.3 

.070 

.104 

.486 

280.0 

235.9 

.110  . 

.096 

.486 

300.0 

239.2 

.160 

.090 

.486 

320.0 

242.5 

.220 

.085 

.486 

340.0 

245.3 

.300 

.081 

.486 

360.0 

247.8 

.390 

.077 

.486 

380.0 

250.2 

.480 

.074 

.486 

400.0 

252.8 

.600 

.072 

.486 

420.0 

255.4 

.730 

.071 

.486 

440.0 

257.9 

.870 

.069 

.486 

460.0 

260.2 

1.030 

.068 

.486 

480.0 

262.6 

1.200 

.067 

.486 

500.0 

264.7 

1.400 

.066 

.486 

520.0 

266.5 

1.630 

.065 

.486 

540.0 

268.3 

1.870 

.064 

.486 

560.0 

270.1 

2.110 

.063 

.486 

580.0 

271.9 

2.300 

.062 

.486 

600.0 

273.8 

2.490 

.061 

.486 

620.0 

275.7 

2.680 

.060 

.486 

640.0 

277.6 

3.290 

.059 

.486 

660.0 

279.3 

4.680 

.059 

.486 

680.0 

281.0 

6.070 

.058 

.486 

700.0 

282.7 

7.460 

.058 

.486 

720.0 

284.2 

8.560 

.058 

.486 

740.0 

285.1 

8.810 

.057 

.486 

760.0 

286.0 

9.050 

.057 

.486 

780.0 

286.8 

9.290 

.057 

.486 

800.0 

287.8 

9.540 

.056 

.486 

820.0 

288.9 

9.990 

.056 

.486 

840.0 

290.1 

10.520 

.055 

.486 

880.0 

291.3 

11.570 

.054 

.486 

900.0 

293.7 

12.090 

.053 

.486 

920.0 

294.8 

12.800 

.052 

.486 

940.0 

295.9 

13.550 

.051 

.486 

Table  3.  (Page  1) 
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PRESSURE  TEMPERATURE 

mb  °K 


W  H20  W  03  W  C02 

G/KG  OG/G  G/ KG 


I 

I 


1013.0 


297.0 

298.1 

299.2 
300.0 


14.300  .050  .486 
15.060  .050  .486 
15.810  .049  .486 

16.300  .048  .486 


Table  3.  Tropical  atmospheric  profile  variables. 


PRESSURE 

mb 

TEMPERATURE 

#K 

W  H20 
G/KG 

W  03 
UG/G 

W  C02 
G/KG 

1.8 

270.0 

.003 

8.600 

.486 

20.0 

229.8 

.003 

9.775 

.486 

40.0 

222.6 

.003 

6.597 

.486 

60.0 

219.0 

.003 

4.260 

.486 

80.0 

217.1 

.003 

2.782 

.486 

100.0 

216.0 

.003 

1.813 

.486 

120.0 

216.0 

.003 

1.258 

.486 

140.0 

216.0 

.003 

.978 

.486 

160.0 

216.0 

.004 

.782 

.486 

180.0 

216.2 

.006 

.675 

.486 

200.0 

220.2 

.010 

.568 

.486 

220.0 

224.3 

.014 

.465 

.486 

240.0 

228.4 

.019 

.365 

.486 

260.0 

231.7 

.030 

.292 

.486 

280.0 

234.8 

.042 

.223 

.486 

300.0 

238.1 

.064 

.189 

.486 

320.0 

241.3 

.086 

.157 

.486 

340.0 

244.0 

.133 

.133 

.486 

360.0 

246.5 

.188 

.112 

.486 

380.0 

249.0 

.239 

.096 

.486 

400.0 

251.6 

.287 

.089 

.486 

420.0 

254.2 

.336 

.082 

.486 

440.0 

256.4 

.396 

.077 

.486 

460.0 

258.3 

.462 

.073 

.486 

480.0 

260.3 

.527 

.069 

.486 

500.0 

262.2 

.599 

.067 

.486 

520.0 

264.0 

.673 

.065 

.486 

540.0 

265.7 

.748 

.063 

.486 

560.0 

267.5 

.841 

.062 

.486 

580.0 

269.1 

.976 

.060 

.486 

600.0 

270.7 

1.111 

.058 

.486 

620.0 

272.4 

1.246 

.057 

.486 

640.0 

273.9 

1.373 

.065 

.486 

660.0 

275.3 

1.495 

.056 

.486 

680.0 

276.8 

1.617 

.055 

.486 

700.0 

278.3 

1.739 

.055 

.486 

720.0 

279.7 

1.909 

.055 

.486 

740.0 

281.0 

2.126 

.055 

.486 

760.0 

282.3 

2.344 

.055 

.486 

780.0 

283.6 

2.561 

.054 

.486 

800.0 

284.9 

2.778 

.054 

.486 

820.0 

285.9 

2.962 

.054 

.486 

860.0 

287.9 

3.322 

.054 

.486 

880.0 

288.9 

3.502 

.054 

.486 

Table  4.  (Page  1.) 


32 


PRESSURE 

mb 

TEMPERATURE 

°K 

W  H20 
G/KG 

W  03 
UG/G 

W  C02 
G/KG 

900.0 

289.9 

3.682 

.054 

.486 

920.0 

290.6 

3.861 

.054 

.486 

940.0 

291.4 

4.039 

.054 

.486 

960.0 

292.1 

4.218 

.054 

.486 

980.0 

292.8 

4.396 

.054 

.486 

1000.0 

293.5 

4.575 

.054 

.486 

1013.0 

294.0 

4.700 

.054 

.486 

Table  4.  Mid-latitude  summer  atmospheric  profile 
variables. 
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1  -  exp(-K  LWC  Az) 


£ CLOUD 

1  “  tcloud  Az 

RT  ,  1 

Z2~  zi  T  pT 

Cloud  Top 
Pressure  (mb) 

Assumed  average  cloud 
ice  or  liquid  water 

(LWC)  content  (gm  ) 

Cloud  penetration  distance 
for  satellite  11  ym  radi¬ 
ance  measurements  (meters) 

100 

0.01 

1660 

200 

0.02 

830 

300 

0.05 

332 

400 

0.10 

166 

500 

0.20 

83 

600 

0.33 

50 

700 

0.50 

33 

800 

1.00 

20 

900 

1.00 

20 

1000 

1.00 

20 

A  mass  absorption  coefficient  (K)  of  0.045 

2  -1 

m  g  was  used  at  all  levels 

thereby  allowing  cloud  emissivity  to  be  determined  by  variations  in 
cloud  water  content  and  cloud  thickness.  (After  Cox  and  Griffith, 
1978) . 
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CLOUD 

TOP 

CLOUD 

BASE 

LUC 

/  3 
g/cm 

COMP. 

T 

COMP. 

a,N«l 

CLOUD 

TOP 

CLOUD 

BASE 

LUC 

,  3 
g/cm 

COMP. 

T 

COMP. 

a,N«l 

100 

200 

.010 

.15 

.85 

560 

660 

.248 

.00 

1.0 

120 

220 

.011 

.15 

.85 

580 

680 

.279 

.00 

1.0 

140 

240 

.013 

.13 

.87 

600 

700 

.313 

.00 

1.0 

160 

260 

.016 

.10 

.90 

620 

720 

.352 

.00 

1.0 

200 

300 

.021 

.08 

.92 

640 

740 

.395 

.00 

1.0 

220 

320 

.025 

.06 

.94 

660 

760 

.443 

.00 

1.0 

240 

340 

.029 

.04 

.96 

680 

780 

.498 

.00 

1.0 

260 

360 

.034 

.03 

.97 

700 

800 

.559 

.00 

1.0 

280 

380 

.040 

.02 

.98 

720 

820 

.628 

.00 

1.0 

300 

400 

.047 

.01 

.99 

740 

840 

.705 

.00 

1.0 

320 

420 

.055 

.01 

.99 

760 

860 

.792 

.00 

1.0 

340 

440 

.064 

.00 

1.0 

780 

880 

.890 

.00 

1.0 

360 

460 

.075 

.00 

1.0 

800 

900 

1.00 

.00 

1.0 

380 

480 

.087 

.00 

1.0 

820 

920 

1.00 

.00 

1.0 

400 

500 

.100 

.00 

1.0 

840 

940 

1.00 

.00 

1.0 

420 

520 

.110 

.00 

1.0 

860 

960 

1.00 

.00 

1.0 

440 

540 

.124 

.00 

1.0 

880 

980 

1.00 

.00 

1.0 

460 

560 

.139 

.00 

1.0 

900 

1000 

1.00 

.00 

1.0 

480 

580 

.156 

.00 

1.0 

920 

1000 

1.00 

.00 

1.0 

500 

600 

.175 

.00 

1.0 

940 

1000 

1.00 

.00 

1.0 

520 

620 

.197 

.00 

1.0 

960 

1000 

1.00 

.00 

1.0 

540 

640 

.221 

.00 

1.0 

980 

1000 

1.00 

.00 

1.0 

Table  5.  Specified  Cloud  Parameters. 
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linear  form  It  is  In  y  ■  In  b  +  mx,  where  m  ■  slope  and  In  b  Is  the  y 
intercept  value.  In  this  model  the  known  independent  variable  x  is 
given  below: 


x  -  L  (p  ,,  a)  =  ct  L  ,,  +  (1-a)  L  , 

V  *w£  veld  ves 

where  is  the  satellite  measured  spectral  radiance  for  a  given  spot, 
a  is  the  effective  cloud  amount,  L  ^  i-8  the  radiance  from  the  cloudy 
area  of  the  spot,  and  L^cs  is  the  radiance  from  the  clear  sky  portion 
of  the  spot.  The  dependent  variable  y  is  equal  to  Pwf  Using  the  data 
represented  in  Figures  6  and  7a  one  can  solve  for  the  equation  para¬ 
meters  m  and  b  for  a  =  1  for  each  wavenumber  line. 

One  could  solve  for  m  and  b  values  for  a  number  of  a  value  lines 
(Figures  6  and  7a  represent  an  a  =  1  line)  using  the  same  procedure. 
However,  an  easier  and  faster  approach  is  to  rewrite  the  linear  regres¬ 
sion  model  in  terms  of  a.  First  note  Figure  7b  which  shows  that  re¬ 
gardless  of  the  a  value  line  plotted,  as  approaches  Lvcg,  pw^ 
approaches  the  constant  value  D.  Furthermore,  since  x  is  a  linear 
combination  using  LVcs  and  a,  one  might  expect  to  be  able  to  write  the 
linear  model  in  terms  of  L,,  and  a  which  in  effect  gives  an  equation 
that  represents  a  family  of  curves  in  pw^  and  a.  In  other  words  con¬ 
sider  the  model 


y  ■  b  exp  (Cx  ■?  a) . 


Analysis  showed  that  for  such  a  model 
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b  -  D  t  exp  (CL  *  a). 

Vcs 

This  results  in  the  following  relationship 

C(L  -  L  ) 

y  =  d  exp [  "  Pwf-  (7) 

In  this  equation  pw£  is  simply  related  to  both  and  a.  Ly  is  meas¬ 
ured  and  L^cs  is  known.  For  the  ot  =*  1  case  plotted  in  Figures  6  and 
7a,  constant  C  is  simply  the  lines'  slope  m,  and  constant  D  =  exp 

(CL  +  In  b)  . 
vcs 

Table  6  gives  the  respective  values  of  C  and  D  and  Lycs  for  each 
wavenumber  band.  Computationally  these  values  are  used  in  Eq.  (7), 
which  is  simply  another  form  of  Eq.  (6)  given  in  Chapter  III.  For  each 
wavenumber  there  is  a  separate  equation  with  two  unknowns  pw£  and  a. 

One  may  solve  a  set  of  two  equations  with  two  unknowns.  Two  wavenumber 
equations  are  chosen.  A  simple  iterative  scheme  that  changes  values  of 
a  in  the  two  equations  is  used  to  solve  simultaneously  for  the  values 
of  p^^  and  a  of  Eq.  (7).  Computationally  one  looks  for  the  best  agree¬ 
ment  between  the  two  calculated  values  of  pw^  as  one  iterates  through 
the  a  values . 

Each  wavenumber  family  of  curves  (Eq.  7)  represents  information 
from  a  part  of  the  atmosphere  as  given  in  Table  6.  As  stated  in  Table 
2  and  illustrated  in  Figures  9  and  10  for  six  spectral  bands  of  inter¬ 
est,  clear  sky  weighting  functions  peak  at  certain  levels  in  the  atmos¬ 
phere.  Most  of  the  radiance  information  at  a  given  wavelength  comes 
from  the  part  of  the  atmosphere  above  this  peak.  A  relatively  small 
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WAVE  NUMBER 

-1 

cm 

WEIGHTING 

LEVEL  ! 

FUNCTION  PEAK 

RANGE  (mb) 

C* 

VARIABLES 

D**  CCR*** 

MID  LATITUDE  SUMMER  ATMOSPHERE 

832.5 

500  - 

690 

3.09 

1007 

.587 

832.5 

690  - 

890 

2.87 

971 

.587 

832.5 

890  - 

995 

3.85 

998 

.587 

747.5 

320  - 

790 

4.55 

920 

.521 

747.5 

790  - 

890 

5.95 

963 

.521 

747.5 

890  - 

995 

9.11 

996 

.521 

727.5 

295  - 

600 

6.47 

815 

.433 

707.5 

190  - 

455 

14.0 

581 

.311 

697.5 

190  - 

295 

31.2 

396 

.259 

TROPICAL  ATMOSPHERE 

832.5 

370  - 

630 

3.21 

942 

.578 

832.5 

630  - 

730 

2.62 

856 

.578 

832.5 

730  - 

935 

4.47 

942 

.578 

747.5 

320  - 

730 

4.12 

795 

.513 

747.5 

730  - 

935 

10.9 

899 

.513 

727.5 

295  - 

660 

6.00 

744 

.428 

707.5 

100  - 

400 

13.9 

600 

.302 

697.5 

100  - 

295 

26.2 

381 

.234 

*  units  of 

[W/m2  str  5 

-1,-1 
cm  J 

**  modified  Intercept  in  units  of  mb. 

explained 

1  in  text 

2 

***  clear  column  radiance  (W/m  str  5 

-1. 
cm  ) 

Table  6.  Cloud  retrieval  empirical  equation  variables. 
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(quu)  ajnss&y 


Figure  9.  Tropical  clear  sky  spectral  weighting  curves 


(qui)  ajnssaJd 
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portion  of  the  total  atmospheric  spectral  radiance  received  at  satel¬ 
lite  level  originates  below  the  clear  sky  weighting  function  peak  level. 
Thus,  any  effective  radiating  surface  below  this  clear  sky  peak  level 
will  not  be  strongly  indicated  in  the  satellite-received  radiance  value. 

This  effect  can  be  seen  in  Figures  6  and  7a.  Above  the  respective 
wave  band  natural  clear  sky  weighting  function  peak,  the  relationship 
of  In  pwf  to  (satellite  received  spectral  radiance)  is  nearly  linear 
with  a  flat  slope.  However,  below  this  point  the  slope  steepens  in¬ 
dicating  very  little  change  in  for  a  change  in  pw^  (level  of  the 
effective  radiating  surface  since  ot  ■  1)  . 

This  effect  is  also  evident  in  Figures  11  and  12  which  plot  wave¬ 
length  dependent  signal  to  noise  ratio  as  a  function  of  height  for  the 
tropical  and  mid- latitude  atmospheres  respectively.  The  signal  - 
LVcs  ^SeC  Eq'  ^  and  t*ie  error  values  are  given  in  Chapter  V.  In  the 
figures,  the  solid  lines  indicate  a  *  1  (overcast  case)  and  the  dashed 
lines  are  for  a  less  than  1.  For  a  values  less  than  1,  the  signal  to 
noise  ratio  is  reduced  by  a  factor  of  a.  For  example,  if  a  =  .5  the 
signal  to  noise  ratio  will  be  half  the  value  it  is  for  a  ■  1.  Note 
that  as  a  ■+■  0  and  as  pw^  (with  constant  o  ■  1)  1000  mb  the  signal  to 

noise  ratio  decreases.  In  other  words  less  information  is  available 
from  which  to  make  curve  shape  retrievals  as  these  boundaries  are  ap¬ 
proached.  However,  as  illustrated  in  Figures  4  and  5  by  the  small 
curve  shape  change  the  VIRES  of  the  atmosphere  is  least  sensitive  to 
errors  near  these  boundaries .  Figures  11  and  12  seem  to  indicate  that 
the  wavenunber  band  curves  peaking  near  the  surface  should  give  the 
best  results  overall.  However,  the  *■  £^2  assumption  limits  the  use 


of  these  widely  spaced  wavelength  curves  higher  in  the  atmosphere  where 
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Figure  11.  Signal  to  noise  ratio  vs.  overcast  effective  cloud  top 
pressure  for  tropical  model.  NOTE:  Signal  =>  LVcg  -  Lv. 


EFFECTIVE  CLOUD  TOP  PRESSURE  (mb) 
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Figure  12.  Signal  Co  noise  ratio  vs.  overcast  effective  cloud  top 
pressure  for  the  mid-latitude  summer  model. 


i 
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nonblack  clouds  are  common.  Also,  the  slope  and  linearity  of  the 
curves  in  Figures  6  and  7a, b  for  these  wavebands  limit  their  use  at 
upper  tropospheric  levels . 

Therefore,  since  the  curve  shape  parameter  radiance  signal  is  re¬ 
lated  to  the  clear  sky  weighting  function  peak  level,  only  the  areas  of 
the  atmosphere  indicated  by  solid  lines  in  Figures  6  and  7a, b  for  the 
respective  wavebands  are  used  in  the  retrieval  calculations.  This  has 
the  effect  of  maximizing  the  useful  shape  parameter  information  con¬ 
tained  in  the  satellite  radiance  values. 

Two  further  points  can  be  made.  First,  the  empirical  RTE  family 
of  curves  approach  requires  linearity  in  In  p  which  restricts  the  use¬ 
ful  part  of  the  curves  in  Figures  6  and  7a, b  as  indicated  in  Table  6. 

To  obtain  best  results  a  few  of  the  solid  lines  plotted  in  Figures  6 
and  7a, b  are  actually  broken  down  into  two  or  three  straight  line  seg¬ 
ments.  All  lines  calculated  resulted  in  the  coefficient  of  determina- 
2 

tion  (r  )  exceeding  0.98.  Secondly,  for  any  given  level  of  the  atmos¬ 
phere  (100  -  980  mb)  there  must  be  at  least  two  different  waveband 
family  of  curves  equations  that  can  be  compared  to  each  other  in  order 
to  determine  the  proper  shape  parameters  pw^  and  a.  Table  6  indicates 
that  this  condition  is  met  and  that  at  times  three  curves  overlap  for 
comparison. 

Neglecting  input  errors  for  now,  it  is  useful  to  evaluate  the 
statistical  characteristics  of  the  previously  described  empirical  re¬ 
trieval.  Figures  13  and  14  show  an  evaluation  of  the  bias  and  SMS 
error  associated  with  retrieval  of  the  two  curve  shape  parameters  pw£ 
and  a.  RMS  error  is  the  standard  deviation  of  retrieval  error.  Bias 
error  is  the  mean  retrieved  value  minus  the  actual  simulated  value. 


Bias  a  RMS  a  Bias  a  RMS 


Effective  Cloud  Top  Categories 


Figure  14.  'No  error'  empirical  model  statistical  analysis  of  o 


Each  point  plotted  is  an  evaluation  of  50  values.  The  50  values  came 
from  the  5  effective  cloud  top  levels  in  a  given  100's  category  (i.e. 
300,  320,  340,  360  and  380  mb)  evaluated  at  10  clear  versus  cloudy 
fractional  weightings  (i.e.  1.0,  0.9,  ...,  0.1).  Values  were  calcu¬ 
lated  for  the  tropical  atmosphere  between  100  and  880  mb  and  between 
200  and  980  mb  for  the  mid-latitude  summer  atmosphere.  Values  above 
the  tropopause  were  not  considered  because  isothermal  conditions  give 

no  indication  of  p^^  changes  in  the  satellite  radiance  values  (i.e. 

3T 

f  0) .  Values  at  900  mb  and  below  in  the  tropical  atmosphere  were 
not  considered  because  the  high  concentration  of  water  vapor  at  these 
levels  essentially  has  the  same  effect  on  the  atmosphere  VIRES  whether 
an  effective  cloud  is  there  or  not. 

From  Figures  13  and  14  it  can  be  seen  that  the  RMS  curves  are  in¬ 
fluenced  by  the  signal  to  noise  ratio  values  discussed  previously.  As 
pointed  out  before,  the  errors  occurring  in  the  retrieval  low  in  the 
atmosphere  have  the  least  effect  on  the  shape  of  the  broadband  weight¬ 
ing  function  curve.  However,  the  error  shown  by  the  bias  curves  is 
only  due  to  the  failure  of  the  empirical  RTE  to  exactly  duplicate  cal¬ 
culation  of  the  RTE  for  every  combination  of  pc  and  a.  The  dotted  line 
shows  an  adjustment  made  to  the  empirical  method  to  remove  a  small  con¬ 
sistently  high  bias.  Variable  D  (an  intercept  value)  of  Eq.  (7)  is 
slightly  reduced  to  nearly  eliminate  this  bias.  Apparently  solving  for 
C  and  D  by  the  method  described  earlier  using  a  =  1  which  is  at  a  bound¬ 
ary,  results  in  this  =  5  mb  bias. 

As  stated  previously  and  as  outlined  in  Figure  la  the  shape  param¬ 
eter  pw£  and  a  are  used  in  the  broadband  computer  routine  referred  to 
earlier  in  this  chapter.  The  program  is  run  using  the  same  atmospheric 
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profile  used  to  compute  the  C  and  D  values  of  Eq.  (7)  and  listed  in 
Table  6.  One  may  run  the  program  with  many  pairs  of  shape  parameters 
to  produce  a  composite  VIRES  curve.  Chapter  VI  gives  examples  of  this 
procedure  and  explains  how  the  VIRES  broadband  curves  are  used  to  yield 
the  proposed  VIRES  index.  The  next  chapter  evaluates  the  retrieval 
errors  under  realistic  system  error  conditions.  Note  also  that  at  the 
end  of  Chapter  V  it  is  shown  that  one  may  produce  useful  composite 
curves  by  combining  the  simulated  spectral  radiances  before  curve  shape 
parameter  retrieval  is  performed. 
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V.  SENSITIVITY  ANALYSIS 

To  help  assess  the  feasibility  of  using  the  atmospheric  VIRES  as  a 
climate  index  in  the  form  of  broadband  infrared  weighting  functions,  it 
is  useful  to  know  the  characteristic  errors  of  inferring  these  profiles 
from  satellite  spectral  radiances.  First,  the  uncertainties  both  ran¬ 
dom  and  bias  associated  with  the  empirical  RTE  technique  used  to  re¬ 
trieve  the  weighting  curve  shape  parameters  pw^  and  a  are  examined. 
Secondly,  the  sensitivity  of  the  final  weighting  curve  product  to  rea¬ 
sonable  errors  is  examined.  All  error  sources  evaluated  are  associated 
with  implicit  or  explicit  assumptions. 

A.  Shape  parameter  retrieval  errors 

Error  sources  examined  in  this  section  include  system- instrument 
errors,  errors  associated  with  temperature,  humidity,  and  profile 
unknowns,  and  errors  due  to  specified  cloud  properties  assumptions. 

The  sensitivity  of  the  weighting  curve  shape  parameters  to  these  likely 
error  sources  is  investigated.  The  chosen  approach  and  specific  error 
magnitudes  are  discussed  below. 

The  calculation  of  the  statistics  used  in  this  analysis  is 
straightforward.  The  pw^  and  a  RMS  errors  are  variations  of  the  re¬ 
trieved  values  about  the  retrieval  mean  pw^  and  a  values.  Thus  the 
RMS  is  a  measure  of  the  retrieved  values  standard  deviation  about  the 
mean.  However,  because  the  RMS  is  large  compared  to  the  mean,  the 
described  RMS  value  is  nearly  equal  to  the  RMS  calculated  the  standard 
way.  Bias  error  is  the  difference  between  the  actual  and  retrieved 


mean  value. 


Important  to  the  calculation  of  the  statistics  is  the  enyirical- 

RTE  technique  boundary  conditions.  Naturally  the  values  of  a  are  found 

between  1.0  and  0.  For  clear  sky  conditions  a  is  given  the  value  0  and 

p  .  *  1000  mb.  Clear  sky  conditions  are  assigned  when  L  >  L  -  2 o  , 
rwf  v  —  vcs  v 

where  v  *  747.5  and  832.5  cm  \  and  is  the  standard  deviation  of  in¬ 
strument  error  at  a  given  wavenumber.  Retrieval  calculations  giving 
radiating  surfaces  above  the  tropopause  are  taken  to  indicate  the 
radiating  surface  location  at  the  tropopause  level  (100  mb  for  tropical 
atmosphere,  and  180  mb  for  mid-latitude  atmosphere).  Because  of  the 
strong  influence  of  water  vapor  in  the  tropical  model,  radiating  sur¬ 
faces  calculated  to  be  below  935  mb  are  taken  to  indicate  clear  sky 
conditions.  Specifics  of  the  error  analysis  are  given  in  the  Appendix. 
A  summary  of  these  results  including  error  assumptions  follows. 

Instrument  system  noise  is  assumed  to  be  Gaussian  with  a  mean  of 

-2  -1 

zero  and  a  standard  deviation  of  0.22  mW  m  sr  cm  for  the  C02  channels 
-2  -1 

and  0.11  mW  m  sr  cm  for  the  window  channel.  These  values  represent 
the  state  of  the  art  precision  of  the  HIRS  instrument  on  TIROS-N 
(Schwalb,  1978).  It  is  also  assumed  that  channel  errors  for  a  specific 
spot  measurement  are  correlated.  Thus  each  simulated  spectral  radi¬ 
ance  value  was  modified  in  a  like  manner  by  an  error  value  distributed 
as  described  above  and  chosen  randomly  for  each  set  of  radiance  values. 
The  effect  of  this  approach,  compared  to  random  assignment  of  individ¬ 
ual  errors  to  the  spectral  radiances  in  a  set,  is  to  give  ~  30%  smaller 
overall  bias  errors  and  s  20%  smaller  RMS  errors  to  the  retrieved  curve 
shape  parameters.  The  bias  and  RMS  errors  are  smaller  because  the  re¬ 
trieval  routine  is  better  able  to  find  a  unique  solution  when  errors 
are  in  the  same  direction  and  of  similar  magnitude.  Real 
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characteristic  sensor  noise  of  course,  lies  somewhere  between  the  two 
descriptions  given  above.  Although  example  comparisons  of  the  results 
from  the  two  approaches  showed  the  trend  described  above,  the  noise 
correlation  approach  was  chosen  for  the  overall  sensitivity  analysis. 

Temperature  and  water  vapor  profile  errors  are  also  assumed  to  be 
Gaussian  with  zero  mean.  Errors  vary  randomly  from  one  20  mb  layer  to 
another  and  temperature  errors  are  assumed  not  to  be  correlated  with 
humidity  errors.  Reasonable  temperature  and  humidity  errors  are  as¬ 
sumed  to  be  £  5°C  RMS  for  temperature  and  an  RMS  £  100%  of  the  correct 
mixing  ratio  (Wielicki  and  Coakley,  1980).  Although  analysis  was  done 
for  a  number  of  reasonable  values,  the  results  reported  on  in  the  Ap¬ 
pendix  are  for  a  2°C  RMS  and  an  RMS  water  vapor  noise  of  50%  (Susskind 
and  Rosenfield,  1980;  Weinreb  and  Crosby,  1977). 

CO ^  profile  maximum  errors  are  assumed  to  be  of  magnitude  1%. 

2 

This  represents  a  value  that  is  of  the  seasonal  change  in  C02  con¬ 
centration  (Bolin  and  Bischof,  1970;  Miller,  1978).  CO^  is  assumed  to 
be  well  mixed  in  the  troposphere  and,  therefore,  the  1%  error  is  ap¬ 
plied  equally  from  layer  to  layer. 

For  all  profile  errors  a  simple  correction  technique  is  sometimes 
useful.  Using  additional  radiative  information  in  the  form  of  a  clear 
column  radiance  correction  (CCRC)  decreases  error  values.  Application 
of  CCRC  is  explained  in  the  Appendix. 

Cloud  radiative  properties  must  be  specified  or  assumed  to  be  of 
a  certain  nature.  For  emitted  radiation  the  most  important  property 
(and  only  one  examined  here)  is  cloud  effective  emittance.  As  Platt 
and  Stephens  (1980)  have  recently  pointed  out,  there  are  scattering  and 
reflection  components  of  effective  cloud  emittance,  especially  for  ice 


clouds.  However,  for  the  cloud  model  used  in  this  research,  two  prop¬ 
erties  determine  the  cloud  effective  emittance;  cloud  thickness  and 
cloud  liquid  water  content  (LWC) .  Two  errors  that  are  most  likely  to 
occur  are  examined.  First  a  100  mb  thick  cloud  is  specified  when  in 
reality  a  20  mb  thin  cloud  is  present.  In  other  words,  cloud  emittance 
is  specified  too  large.  Second,  a  cloud  of  higher  LWC  than  specified 
is  present.  This  is  the  case  of  cumulonimbus  cloud  at  cirrus  level. 

In  this  situation  cloud  emittance  is  specified  too  small.  This  occurs 
only  above  300  mb  since  below  thin  level  cloud  emittance  is  specified 
at  the  maximum  value  of  1.  As  stated  previously  and  as  shown  by  this 
analysis,  the  retrieval  process  is  essentially  independent  of  cloud 
radiative  properties  specified. 

Table  7  gives  a  summary  of  the  error  sources  and  the  corresponding 
bias  and  RMS  order  of  magnitude  error  for  the  two  curve  shape  param¬ 
eters  retrieved.  The  combination  error  source  includes  a  random  2°C 
RMS  and  50%  RMS  water  mixing  ratio  error.  It  also  includes  instrument 
noise  as  described  above.  Cloud  depth  is  allowed  to  randomly  vary  from 
20  to  180  mb  and  cloud  LWC  varies  randomly  from  +  50%  of  the  specified 
value.  This  table  shows  results  that  testify  to  the  capability  of  this 
retrieval  method. 

B.  Error  effects  on  the  weighting  curves 

The  next  source  of  error  evaluated  is  the  assumption  that  the 
satellite  sensor  views  only  one  effective  cloud  layer  in  non-clear 
situations.  Direct  sensitivity  of  the  weighting  curve  to  this  error  is 
investigated.  The  error  is  examined  by  means  of  a  sensor  resolution 
argument.  We  compare  two  models.  In  the  first  model  the  sensor 
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Error 

Sources* 

Pwf  bias  (mb) 

Statistical  Values* 
Pw£  RMS  (mb)  ct  bias 

O  RMS 

Empirical  RTE 

0 

5 

±  .01 

.03 

Sensor  noise 

±  15 

80 

±  .04 

.15 

Temperature  +  humidity 

±  30 

160 

±  .06 

.30 

Temperature  +  humidity 
with  CCRC* 

±  15 

20 

±  .06 

.05 

C02  profile 

±  30 

40 

t  .02 

.02 

C02  profile  with  CCRC 

±  2 

10 

±  .01 

.02 

Specified  emittance 

±  2 

8 

±  .01** 

.03** 

Coni)  inat  ion 

±  15 

160 

±  .05 

.35 

Combination  with  CCRC 

±  15 

90 

±  .05 

.20 

*  defined  in  text 

**  not  plotted  In  the  Figures 


Table  7.  List  of  representative  errors  in  shape  parameter  retrieval 
due  to  given  error  sources. 
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detects,  for  example,  200  sets  of  radiances,  representing  single  layer 
effective  clouds.  This  model  is  compared  to  the  second  model  which 
senses  100  sets  of  radiance  values  for  the  same  area  of  effective  cloud. 
The  second  set  of  radiances  are  simply  averages  of  the  first  set  of 
radiances  taken  by  pairs.  For  both  cases  curve  shape  parameters  are 
retrieved  assuming  single  layer  effective  clouds  and  composite  weight¬ 
ing  function  curves  which  represent  the  same  200  spot  (high  resolution) 
area  are  computed.  The  two  curves  are  statistically  compared  to  test 
the  one  layer  error  assumption.  In  the  200  spot  specified  model  the 
assumption  is  totally  true,  while  in  the  100  pair  model  the  assumption 
is  completely  false. 

In  this  analysis  the  200  spot  area  represents  a  sample  from  an 
effective  cloud  layer  distribution.  Three  different  distributions  of 
effective  cloud  top  heights  are  specified  to  provide  a  more  complete 
test.  The  test  in  general  as  well  as  the  distributions  are  described 
in  detail  in  the  Appendix.  The  statistical  test  used  is  described  in 
Chapter  VI. 

Table  8  shows  the  results  of  the  single  effective  cloud  layer 
assumption  analysis.  The  statistical  test  values  indicate  that  for  the 
situations  of  curve  compositing  the  errors  associated  with  the  single 
layer  assumption  are  minimized.  For  the  test  cases  presented,  the  re¬ 
trieved  composite  weighting  function  curves  are  extraordinarily  similar 
to  the  true  (specified)  curves.  Although  the  statistical  test  is  de¬ 
signed  to  indicate  curve  differences,  small  T  values  imply  likeness. 

In  other  words,  if  the  calculated  T  values  are  as  large  or  larger  than 
the  critical  T  value  then  HQ  (the  curves  are  the  same)  can  be  rejected 
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Table  8.  Statistical  test  of  the  differences  between  the  curves  in  Figures  40-42. 


at  the  given  P  value  significance  level.  Otherwise,  is  accepted  as 
true,  the  differences  in  the  curves  being  explained  by  chance  happen¬ 
ings. 

This  last  section  looks  at  the  sensitivity  of  the  final  product 
broadband  weighting  curve  to  the  factors  that  are  most  important  to 
curve  parameter  retrieval,  temperature  and  humidity.  The  two  climato¬ 
logical  profiles  (Tables  3  and  4)  are  assumed  to  be  correct  when  using 
the  retrieval  method,  i.e.  when  determining  C  and  D  values  and  when 
calculating  the  broadband  curves  from  the  retrieved  curve  shape  param¬ 
eters.  However,  each  of  these  profiles  are  modified  to  produce  'actual 
cases .  Temperature  and  humidity  at  each  level  in  the  climatological 
profiles  are  increased  or  decreased  by  a  specified  percentage  to  pro¬ 
duce  these  'actual'  cases.  From  these  'actual'  profiles  the  simulated 
satellite  radiances  are  calculated.  Finally,  the  broadband  weighting 
curve  retrieved  assuming  a  climatological  atmosphere  (without  using 
CCRC)  when  an  'actual'  profile  exists  is  compared  to  the  weighting 
curve  produced  directly  from  the  true  atmospheric  profile  and  'true' 
specified  cloud  distribution.  The  specified  cloud  profile  used  is 
called  1N24  and  it  represents  the  actual  24  hour  average  cloud  condi¬ 
tions  for  the  B  array  north  section  during  Phase  III  of  GATE  (see  Cox 
and  Griffith,  1978).  This  cloud  distribution  as  well  as  those  used 
in  Chapter  VI  and  the  Appendix  are  given  in  Table  9. 

Changing  the  profiles  has  the  following  general  effects.  In¬ 
creasing  the  temperature  or  decreasing  the  humidity  in  the  profiles 
gives  radiance  values  higher  than  those  calculated  using  a  climato¬ 
logical  profile.  Likewise  given  specified  pw^  and  a  values  the  radi¬ 
ances  from  a  climatological  profile  will  be  higher  than  those 
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Eff. 

Cld. 
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Table  9.  Model  effective  cloud  distributions  (frequency  of  occurrences). 
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calculated  using  cooler  or  moister  profiles.  In  general  then,  the  pwf 
values  retrieved  using  radiances  from  a  warmer  or  drier  atmosphere  will 
be  a  few  millibars  (20  -  40  mb  given  the  conditions  in  the  next  para¬ 
graphs)  larger  (lower  in  the  atmosphere)  than  the  specified  value.  Just 
the  opposite  is  true  for  the  cooler  or  moister  case.  Typically,  re¬ 
trieved  a  values  remain  equal  or  sligihtly  smaller  than  the  specified 
a  value. 

To  illustrate  the  findings.  Table  10  gives  the  results  of  this 
analysis  for  the  case  of  increasing  the  temperature  of  the  climato¬ 
logical  profile  by  1%  at  each  level  and  for  the  case  of  increasing  the 
water  mixing  ratio  at  each  level  by  20%.  Results  for  both  tropical  and 
mid- latitude  cases  are  given. 

It  is  clear  that  for  these  cases  the  VIRES  index  (described  in  the 
next  chapter)  is  not  sensitive  to  these  reasonable  assumption  errors. 

The  small  T  statistic  values  also  indicate  that  the  broadband  weighting 
curves  are  little  effected.  Looking  at  the  specified  and  retrieved 
effective  cloud  distributions  one  notices  good  agreement  for  the  high, 
middle  and  combined  low  and  clear  categories  but  poor  agreement  for  the 
low  and  clear  categories  taken  separately.  This  simply  supports  what 
has  already  been  said  concerning  poor  retrieval  capabilities  in  the  low 
troposphere.  It  also  shows  the  advantage  of  retrieving  the  VIRES  index. 

A  qualifying  statement  is  in  order.  A  more  severe  test  of  this 
retrieval  method  would  be  to  not  only  shift  the  temperature  and  humid¬ 
ity  profiles  but  to  also  change  their  shape.  This  was  done  of  course, 
for  the  tests  in  Section  A.  Also,  a  rather  flat  distribution  of  clouds 
with  some  clear  sky  is  a  reasonable  test,  but  poorer  results  would  be 
expected  from  a  sharper  (more  peaked)  distribution.  Different  shaped 
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effective  cloud  distributions  and  their  respective  VIRES  indexes  are 
examined  in  the  next  chapter. 

The  Appendix  has  a  detailed  description  of  the  error  Investigation 
as  well  as  a  discussion  on  minimizing  the  errors.  For  example,  the 
use  of  CCRC  improves  the  accuracy  of  the  retrieval  over  using  pre¬ 
selected  climatological  profiles  (Crutcher  and  Mesere,  1970;  Jenne 
et  al.  1974;  Smith  et  al.  1972).  Nevertheless,  the  true  strength  of 
this  climate  index  approach  is  that  radiation  information  is  used  to 
produce  an  index  directly  from  radiation  theory.  Thus,  the  error 
prone  process  of  inferring  standard  physical  meteorological  variables 
is  avoided. 
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VI.  WEIGHTING  CURVES  AS  A  CLIMATE  INDEX 

A  first  step  in  developing  a  climate  index  from  composite  weight¬ 
ing  curves  is  to  find  an  objective  way  to  infer  from  the  curves  when 
there  are  different  types  of  VIRES  (indicators  of  climatological  cloudi¬ 
ness).  Section  A  of  this  chapter  describes  a  statistical  technique  and 
shows  that  it  is  a  powerful  test  for  objectively  determining  when  two 
weighting  function  curves  and  thus  two  effective  cloud  top  distribu¬ 
tions  are  different.  The  remainder  of  the  chapter  deals  with  developing 
and  using  a  climate  index  based  on  the  retrieved  broadband  weighting 
function  curves. 

A.  Sign  Test  with  Fisher’s  Method 

Utilizing  the  broadband  infrared  RTE  described  earlier,  the  weight¬ 
ing  function  curves  are  computed  using  a  finite  differencing  scheme 
which  provides  one  curve  value  every  20  mb  between  100  and  1000  mb. 

These  45  point  values  describe  the  weighting  function  in  the  part  of 
the  atmosphere  influenced  by  effective  clouds.  These  curves  quite 
obviously  represent  no  standard  statistical  distribution.  Therefore, 
when  comparing  two  broadband  weighting  function  curves  one  must  use 
non-parametric  statistics. 

Because  the  continuous  curve  is  computed  using  discrete  points, 
it  is  easy  to  compare  two  curves  by  pairing  their  respective  point 
values.  A  non-parametric  statistical  test  may  then  be  used  to  test  the 
likeness  or  difference  of  the  paired  values.  Although  the  Wilcoxon 
Matched  Pair  Signed  Rank  Test  and  the  Kolomogorov  -  Smirnov  Test  are 
sometimes  useful  in  such  cases,  the  simple  Sign  Test  gives  the  best 
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results  In  this  case.  The  following  is  a  description  of  the  Sign  Test 
(Beyer,  1971). 

In  this  test,  observations  from  sample  x  and  sample  y  are  paired 
and  the  differences  are  calculated.  The  null  hypothesis  (H^)  is  that 
the  difference  d^  has  a  distribution  with  median  zero,  i.e.  the  true 
proportion  of  positive  (negative)  signs  is  equal  to  P  =  The  proba¬ 
bility  of  positive  (negative)  signs  is  given  by  the  binomial  probabil¬ 
ity  function  f(x)  such  that 

f(x)  -  f(x;  n,  P  =  -|)  -  (“)  (|)  , 

where  x  is  either  the  number  of  positive  signs  or  negative  signs,  and 
n  is  the  sum  of  positive  and  negative  signs.  The  probability  P±  of  the 
samples  being  the  same  (true  R_)  given  x  positive  (negative)  signs  is 
given  by  the  expression 


Pt  (x  <  k)  -  2[Z  (")]  (■ |)“ , 

x=0 

where  P^x)  is  the  Sign  Test  Probability  Function  for  the  two-tailed 
case.  The  test  is  two-tailed  because  there  is  no  reason  to  expect  more 
positive  signs  than  negative  signs.  For  example,  if  k  =  3  and  n  *  15, 
the  probability  of  x  being  three  or  less  (either  of  positive  or  nega¬ 
tive  sign)  if  Hq  is  true  is  given  by  the  following  calculation. 


2  id5)  d5)  d5)  d5)]  d) 


p.  (x  <  3) 


.035156. 
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In  other  words,  1  -  .035156  *  .964844  is  the  probability  that  is 
false.  In  this  case  Hq  would  be  rejected  at  the  5%  level.  When  they 
occur,  zero  differences  are  excluded  because  they  contribute  no  in¬ 
formation  for  deciding  whether  positive  or  negative  differences  are 
more  likely. 

To  evaluate  the  usefulness  of  the  Sign  Test  sixteen  effective 
cloud  top  distributions  are  specified.  Figure  15  and  Table  9  describe 
these  distributions  in  detail.  Note  that  three  of  these  distributions 
are  used  for  the  error  analysis  in  the  last  chapter  and  detailed  in  the 
Appendix.  The  specified  distributions  are  intended  to  represent  the 
type  of  observed  effective  cloud  top  distributions  illustrated  in  Fig¬ 
ure  16.  The  corresponding  infrared  broadband  weighting  function  curves 
for  these  sixteen  distributions  are  given  in  Figures  17-23  for  the 
tropical  atmosphere.  These  are  the  curves  we  wish  to  discriminate  be¬ 
tween. 

The  characteristic  shape  of  the  curves  in  Figures  17-23  result  in 
two  problems  when  using  the  Sign  Test.  First,  the  curves  exhibit  a 
characteristic  known  as  crossover.  When  comparing  two  curves  such  as 
FLATS  and  BIMOS,  or  FLATS  and  1PK3S  in  Figure  17,  one  notices  the 
curves  crossover  each  other  many  times  in  the  first  case  and  at  least 
once  in  the  second  case.  Because  the  curves  indicate  vertical  weight¬ 
ing  and  total  weighting  changes  little  from  case  to  case,  the  curve 
that  shows  large  weighting  at  one  vertical  position  must  show  small 
compensating  weighting  at  another  point  in  the  vertical.  This  con¬ 
straint  is  the  primary  cause  of  curve  crossovers.  The  crossover  effect 
reduces  the  ability  of  the  Sign  Test  to  discriminate  between  two  ob¬ 
viously  different  curves.  This  is  clearly  seen  in  Figure  19.  1PK2 
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Figure  15.  Plots  of  effective  cloud  top  model  distributions 
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Figure  15  (Continued) .  Plots  of  effective  cloud  top  model  distribu- 


66 


tr-trtrt 


k  t  fe 
k  k  k  k  k 


.fc. 

k 

b  b  bk- 

kkfcb 

1  1 
^  ,  c 

»^Ub-b 

I  I  ' 


br  kr  1 


El  &=.  §E  t. 

1  i  1  '  iiii 

pr-  |l_  |£-  F^g=-_  g 

k_  L.  1  !  ’"" 


fe.  fe_  St  hr  a 


t  t  k  i 


i  ! 


fcrfcr 


U.  L_ 

i  i 


Figure  16.  Unpublished  figure  by  G.  G.  Campbell  showing  histograms 
(5°K  resolution,  16  cat.)  of  geosynchronous  satellite  IR 
window  radiance  temperatures  that  approximate  the  cloud 
top  distribution.  Each  histogram  represents  200  x  200  km 
area.  The  total  figure  represents  data  taken  on  16  Nov. 
1978  at  1500  L  for  the  Pacific  Ocean  region  10°N  -  20°S 
latitude  and  235°E  -  265°E  longitude. 

NOTE:  The  length  of  the  x-axis  represents  a  50%  frequency  and  the 
y-axis  has  temperature  decreasing  up  the  axis. 


Broadband  weighting  curves  for  given  specified  distributions 


Broadband  weighting  curves  for  given  specified  distributions 


22.  Broadband  weighting  curves  for  given  specified  distributions 


23.  Broadband  weighting  curves  for  given  specified  distributions. 
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and  1PK6  are  compared.  For  this  case  n  ■  45  and  x  ■  19.  The  Sign  Test 
cannot  reject  Hq  even  at  the  25%  level.  The  test  shows  that  there  is 
greater  than  1  chance  in  4  that  these  two  curves  are  the  same.  The 
crossover  effect  is  to  be  blamed  for  this  test  deficiency. 

The  crossover  effect  may  be  greatly  reduced  by  considering  the 
curves  in  three  equal  parts  of  15  pairs  each  instead  of  taken  as  a 
whole  (45  pairs).  The  curves  may  be  divided  into  three  physically 
distinct  parts.  This  division  parallels  the  well  known  high,  middle 
and  low  cloud  categories .  The  part  between  980  and  700  mb  represents 
the  lower  tropospheric  area  of  high  water  vapor  content  (approximately 
85%  of  the  total) .  In  this  area  water  vapor  represents  an  important 
source  of  infrared  emission  as  illustrated  by  the  low  level  weighting 
function  peak  for  clear  sky  in  Figure  17.  From  380  to  100  mb  another 
naturally  occurring  clear  sky  weighting  function  peak  is  seen.  This  is 
also  the  area  of  cirrus  cloud  formation.  Since  ice  clouds  are  somewhat 
different  emittors  than  water  clouds  it  is  not  unreasonable  to  examine 
their  respective  influence  on  the  weighting  curves  separately.  Final¬ 
ly,  the  680  to  400  mb  layer  is  found  between  the  two  layers  described 
above.  It  lacks  the  consistent  large  amount  of  water  vapor  and  the 
clear  sky  weighting  function  peaks.  It  is  also  true  that  the  emitting 
structure  within  these  three  divisions  are  for  the  most  part  determined 
using  different  pairs  of  spectral  radiative  values. 

Fisher’s  Method  of  Combining  Three  P  Values  is  used  with  the  Sign 
Test  when  comparing  corresponding  pairs  from  two  curves  in  three  parts 
(Fisher,  1958).  Fisher  found  that  the  probability  function  distribu¬ 
tion  of  -  2  times  the  sura  of  the  logarithms  of  m  independent  P  values 
is  Chi-Squared  with  2m  degrees  of  freedom.  Thus,  the  three  P  statistic 


values  are  distributed  as  Chi-Squared  with  6  degrees  of  freedom.  The 
calculation  of  the  T  statistic  is  done  using  the  values  in  Table  11 
computed  from  the  following  equation. 


T  -  (-2)  [In  ?1  +  In  P2  +  In  P3]  =  X2d  m 

For  example,  if  In  P^  ■  -4.9,  In  P2  ■  -6.3,  and  In  P^  =  -2.1  (see  Table 

2 

11),  the  T  statistic  is  calculated  as  26.6.  Since  for  P  «*  .01,  ■ 

16.8,  one  would  reject  Hq  at  the  1%  level  (i.e.  less  than  1  chance  in 
100  that  the  two  samples  are  the  same  and  Hg  is  true) .  The  critical 
T  values  for  a  specified  P  level  are  given  in  Table  11. 

Table  12  gives  the  T  statistic  for  the  different  curve  comparisons. 
By  considering  a  critical  value  of  T  =  10.6  (P  «=  0.10),  Table  12  shows 
that  the  test  indicates  curve  differences  in  all  cases  except  one 
(BIM01  vs.  1PK3S).  Notice  also  that  the  test  says  FLAT1  vs.  FLATS  are 
different  curves  when  in  fact  they  are  nearly  identical.  These  two 
cases  point  up  the  second  drawback  of  this  test.  Although  the  test 
allows  for  zero  differences,  the  43  curve  point  pairs  never  give  an 
absolute  difference  value  of  zero  primarily  because  values  are  computed 
beyond  reasonable  significant  figures.  Therefore,  at  times  their  dif¬ 
ference  is  very  small,  but  as  computed  for  Table  12  there  are  always  3 
sets  of  15  sign  values  calculated  with  no  zeros.  Notice  that  above 
100  mb  (Figures  17-20)  curve  points  are  identical  but  they  are  not  used 
in  the  Sign  Test. 

To  correct  this  deficiency,  when  the  difference  between  curve 
points  is  smaller  than  a  specified  significance  level,  it  is  set  to 
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n  ” 

15 

X 

InP 

0 

-9.70A0605 

1 

-6.99601 

2 

-A. 90827 

3 

-3.3A795 

A 

-2.13310 

5 

-1.198131 

6 

-  .A9883 

7 

0 

n  * 

12 

0 

1 

2 

3 

A 

5 

-7.62A62 

-5.05967 

-3.255171 

-1.92A18 

-  .9A7536 

-  .2556A9 

6 

0 

n  " 

9 

0 

1 

2 

3 

-5.5A518 
-3.2A26 
-1.7165A 
-  .6776A 

A 

0 

n  => 

6 

0 

1 

2 

-3.A657A 
-1.5198 
-  .37A7 

3 

0 

n  =* 

3 

0 

-1.3863 

1 

0 

n  -  1A 


X 

InP 

0 

-9.010913 

1 

-6.302863 

2 

-A. 35695 

3 

-2.8603106 

A 

-1.717896 

5 

-  .858A3 

6 

-  .23521 

7 

0 

n  = 

11 

0 

-6.931A7 

1 

-A.AA66 

2 

-2.72678 

3 

-1.A8A7 

A 

-  .59997 

5 

0 

n  = 

8 

0 

-A. 8520 

1 

-2.65A8 

2 

-1.2411 

3 

-  .319A3 

n  = 

5 

0 

-2.7726 

1 

-  .9808 

2 

0 

n  = 

2 

0 

-  .6931 

0 

n  * 

13 

X 

InP 

0 

-8.31777 

1 

-5.67871 

2 

-3.79598 

3 

-2.38287 

A 

-1.32108 

5 

-  .5A2910A 

6 

0 

n  ■ 

10 

0 

-6.2383 

1 

-3.8A0A3 

2 

-2.2197 

3 

-1.0678A 

A 

-  . 282A8 

5 

0 

n  = 

7 

0 

-A. 1589 

1 

-2.079A 

2 

-  .79159 

3 

0 

n  = 

A 

0 

-2.079AA 

1 

-  . A700 

2 

0 

n  - 

1 

0  0 


Table  11.  The  logarithm  of  Sign  Test  probabilities  for  given  n  and  x 
values. 


77 


DISTRIBUTIONS 

100  - 
Pos. 

380  mb 
Neg. 

SIGN 

400  - 
Pos . 

TEST 

-  680  mb 
Neg. 

700  - 
Pos . 

980  mb 
Neg. 

T-statistic 

FLATS  vs.  BIMOS 

10 

5 

12 

3 

4 

11 

13.4 

FLATS  vs.  BIM01 

10 

5 

15 

0 

8 

7 

21.8 

FLATS  vs.  1PK3S 

6 

9 

10 

5 

15 

0 

22.8 

FLAT1  vs.  FLATS 

6 

9 

12 

3 

1 

14 

21.7  + 

FLAT1  vs.  BIMOS 

10 

5 

12 

3 

4 

11 

13.4 

FLAT1  vs.  BIM01 

10 

5 

15 

0 

8 

7 

21.8 

FLAT1  vs.  1PK2 

4 

11 

15 

0 

15 

0 

43.1 

FLAT1  vs.  1PK3S 

5 

10 

11 

4 

15 

0 

26.1 

FLAT1  vs.  1PK4 

10 

5 

6 

9 

15 

0 

22.8 

FLAT1  vs.  1PK6 

15 

0 

5 

10 

12 

3 

28.5 

FLAT1  vs.  1PK8 

15 

0 

12 

3 

14 

1 

40.1 

BIMOl  vs.  BIMOS 

5 

10 

13 

2 

9 

6 

13.2 

BIM01  vs.  1PK2 

2 

13 

15 

0 

10 

5 

31.6 

BIMOl  vs.  1PK3S 

5 

10 

4 

11 

9 

6 

7.7  + 

BIMOl  vs.  1PK4 

10 

5 

3 

12 

15 

0 

28.5 

BIMOl  vs.  1PK6 

11 

4 

0 

15 

13 

2 

33.5 

BIMOl  vs.  1PK8 

11 

4 

0 

15 

0 

15 

43.1 

1PK2  vs.  1PK3S 

10 

5 

0 

15 

3 

12 

28.5 

1PK2  vs.  1PK4 

11 

4 

0 

15 

14 

1 

37.7 

1PK2  vs.  1PK6 

13 

2 

0 

15 

6 

9 

30.2 

1PK2  vs.  1PK8 

13 

2 

0 

15 

0 

15 

48.6 

1PK3S  vs.  1PK4 

15 

0 

3 

12 

15 

0 

45.5 

1PK3S  vs.  1PK6 

15 

0 

5 

10 

7 

8 

21.8 

1PK3S  vs.  1PK8 

15 

0 

6 

9 

0 

15 

40.1 

1PK4  vs.  1PK6 

•s 

8 

7 

14 

1 

16.4 

1PK4  vs.  1PK8 

5 

9 

6 

0 

15 

22.8 

1PK6  vs.  1PK8 

0 

15 

15 

0 

1 

14 

52.8 

+  referred  to  in  the  text 

Table  12.  Sign  Test  results  without  zero  cutoff  correction. 
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zero.  Figure  18  illustrates  the  small  differences  between  curves  BIM01 
and  1PK3  below  850  mb.  These  differences  are  set  to  zero.  We  refer  to 
this  procedure  as  zero  cutoff  correction.  Table  13  which  includes  more 
curve  comparisons  than  Table  12  gives  the  T  values  with  this  correction 
applied.  Notice  that  most  of  the  change  occurs  in  the  700-980  mb  group. 
Also  notice  that  BIM01  and  1PK3S  now  pass  the  test  for  being  different 
while  FLATS  and  FLAT1  as  well  as  1PK8  and  1PK9  have  small  T  values  that 
indicate  likeness.  It  is  also  instructive  to  note  that  the  test  does 
not  catch  the  difference  between  FLAT1  and  BIMOS  at  the  P  =*  5%  level 
but  it  does  at  the  P  *  10%  level.  As  Figure  15  illustrates,  BIMOS  re¬ 
presents  a  flatter  bimodal  effective  cloud  top  distribution  than  BIM01 
which  explains  why  FLAT1  and  BIM01  are  judged  to  be  different. 

From  the  P  values  given  in  Table  8  it  is  clear  that  T  values 
greater  than  about  10  represent  rejection  of  the  null  hypothesis. 

Curves  with  smaller  T  values  indicate  acceptance  of  H^:  the  curves  are 
the  same.  Referring  to  Table  13,  values  of  T  greater  than  10  indicate 
unlike  curves.  Of  the  48  combinations  of  curve  pairings  only  6  cases 
resulted  in  T  values  where  Hq  cannot  be  rejected  at  the  5%  level  (two 
cases  are  mentioned  above).  In  two  of  these  cases  (BIM02  vs.  BIM03  and 
1PK3  vs.  1PK3M)  the  effective  cloud  top  distributions  are  nearly  the 
same,  explaining  the  statistical  results.  In  the  other  two  cases 
(BIM01  vs.  1PK3  and  1PK8  vs.  1PK7)  the  large  amount  of  water  vapor  in 
the  lower  troposphere  effectively  acts  as  a  radiating  surface  resulting 
in  broadband  weighting  curve  pairs  that  are  quite  similar,  yet  repre¬ 
sent  different  effective  cloud  top  distributions. 

Two  other  points  can  be  made  from  Table  13.  The  test  is  able  to 
distinguish  between  two  slightly  different  distributions  peaking  at 
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DISTRIBUTIONS 

100  - 
Pos. 

380  mb 
Neg. 

SIGN 
400  ■ 
Pos. 

TEST 
-  680  mb 
Neg. 

700  - 
Pos . 

980  mb 
Neg. 

T-statistlc 

FLATS 

vs. 

BIMOS 

9 

5 

10 

2 

0 

5 

13,8 

FLATS 

vs. 

1PK3S 

5 

6 

8 

4 

11 

0 

15.8 

FLATS 

vs. 

BIMOl 

10 

5 

13 

0 

3 

6 

20.4 

FLAT1 

vs. 

FLATS 

5 

3 

0 

0 

0 

0 

0.6  + 

FLAT1 

vs. 

BIMOS 

7 

5 

11 

3 

0 

5 

11.3  + 

FLAT1 

vs. 

1PK3S 

5 

8 

9 

4 

10 

0 

16.2 

FLAT1 

vs. 

1PK2 

3 

11 

15 

0 

11 

0 

39.0 

FLAT1 

vs. 

1PK3 

5 

9 

11 

3 

12 

0 

22.7 

FLAT1 

vs. 

1PK4 

10 

5 

6 

9 

12 

0 

18.6 

FLAT1 

vs. 

1PK5 

13 

0 

1 

12 

12 

0 

43.2 

FLAT1 

vs. 

1PK6 

13 

0 

2 

9 

7 

3 

24.2 

FLAT1 

vs. 

1PK8 

13 

0 

8 

2 

12 

0 

36.3 

FLAT1 

vs. 

BIMOl 

10 

5 

14 

0 

3 

6 

20.4  + 

FLAT1 

vs. 

BIM02 

10 

5 

7 

7 

10 

1 

11.3 

BIMOl 

vs. 

BIMOS 

5 

5 

9 

0 

6 

5 

11.1 

BIMOl 

vs  • 

BIMD2 

10 

5 

0 

10 

7 

0 

23.2 

BIMOl 

vs. 

BIM03 

10 

5 

2 

12 

8 

0 

20.8 

BIM02 

vs. 

BIM03 

0 

0 

3 

6 

4 

0 

5.5  + 

BIMOl 

vs. 

1PK2 

0 

13 

11 

0 

8 

0 

40.2 

BIMOl 

vs. 

1PK3 

5 

5 

6 

6 

8 

0 

9.7  + 

BIMOl 

vs. 

1PK4 

10 

5 

3 

10 

8 

0 

16.9 

BIMOl 

vs. 

1PK5 

10 

4 

1 

14 

8 

0 

27.1 

BIMOl 

vs. 

1PK6 

11 

3 

0 

15 

6 

2 

27.6 

BIMOl 

vs. 

1PK8 

11 

3 

0 

15 

0 

13 

41.8 

1PK2  vs. 

1PK3 

10 

5 

0 

9 

0 

0 

13.5 

1PK2  vs. 

1PK4 

10 

4 

0 

14 

0 

0 

21.5 

1PK2  vs. 

1PK5 

12 

2 

0 

15 

0 

2 

29.5 

1PK2  vs. 

1PK6 

12 

2 

0 

15 

0 

6 

35.1 

1PK2  vs. 

1PK8 

12 

2 

0 

15 

0 

13 

44.8 

Table  13.  (Page  1) 
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DISTRIBUTIONS 

100  - 
Pos. 

380  nb 
Neg. 

SIGN  TEST 

400  -  680  mb 
Pos.  Neg. 

700  - 
Pos. 

980  mb 
Neg. 

T-statistlc 

1PK3  vs .  1PK4 

15 

0 

0 

13 

0 

0 

36.0 

1PK3  vs.  1PK5 

15 

0 

3 

12 

0 

3 

28.9 

1PK3  vs.  1PK6 

15 

0 

4 

11 

0 

7 

32.0 

1PK3  vs.  1PK8 

15 

0 

4 

10 

0 

13 

39.5 

1PK4  vs.  1PK5 

10 

0 

6 

9 

0 

3 

16.2 

1PK4  vs.  1PK6 

10 

0 

8 

6 

0 

7 

21.3 

1PK4  vs.  1PK8 

10 

0 

9 

5 

0 

13 

30.8 

1PK5  vs.  1PK6 

3 

0 

11 

4 

0 

7 

15.4 

1PK5  vs.  1PK8 

3 

0 

13 

2 

0 

13 

29.2 

1PK6  vs.  1PK7 

0 

0 

11 

0 

0 

10 

26.3 

1PK6  vs.  1PK8 

0 

0 

11 

0 

1 

12 

25.2 

1PK8  vs.  1PK7 

0 

0 

0 

5 

8 

5 

6.6  + 

1PK8  vs.  1PK9 

0 

0 

0 

0 

4 

5 

0.0  + 

1PK3  vs.  1PK3S 

5 

0 

2 

11 

0 

4 

17.3  + 

1PK3  vs.  1PK3M 

10 

5 

8 

2 

0 

0 

6.8  + 

1PK3S  vs.  1PK3M 

10 

5 

12 

2 

6 

0 

18.0  + 

BIM01  vs.  1PK3S 

5 

5 

3 

10 

8 

0 

14.5  + 

1PK3  vs.  1/2PK3 

15 

0 

4 

9 

0 

15 

41.5  + 

1PK5  vs.  1/2PK5 

3 

0 

13 

2 

0 

15 

32.0  + 

+  Referred  to  in  text 


Table  13.  Sign  Test  with  zero  cutoff  correction. 
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the  same  level.  For  example,  consider  the  cases  1PK3  vs.  1PK3S  and 
1PK3S  vs.  1PK3M.  Table  9  gives  the  cloud  frequency  distribution  for 
these  two  similar  cases  of  a  single  peak  within  the  300  mb  category. 
The  1PK3  vs.  1PK3S  case  represents  only  a  15%  cloud  amount  difference 
in  the  300  mb  category.  It  is  also  clear  that  the  test  is  most  sensi¬ 
tive  when  total  300  mb  category  cloud  amounts  are  less  than  50%. 
Secondly,  the  test  can  distinguish  between  two  cases  where  the  effec¬ 
tive  cloud  top  distribution  is  the  same  but  the  proportion  of  clear 
versus  cloudy  varies.  For  example,  1PK3  represents  an  overcast  case 
whereas  y  PK3  represents  ^  clear  sky.  The  test  is  also  done  using 
1PK5  and  PK5.  In  summary  these  results  show  chat  the  statistical 
test  is  capable  of  distinguishing  between  VIRES  curves  representing 
either  cloud  distribution  shifts  in  the  vertical  or  in  total  amount. 

It  is  successful  in  about  90%  or  more  of  all  cases  at  the  95%  confi¬ 
dence  level. 

The  zero  cutoff  correction  requires  further  explanation.  The 
small  difference  value  is  considered  zero  when  the  difference  is  less 
than  a  specified  constant.  This  constant  is  5%  of  the  average  maximum 
point  value  of  all  points  between  100  and  980  mb  for  all  the  weighting 
curves  plotted.  In  reference  to  Figures  40-43,  this  is  approximately 
a  distance  of  .02  on  the  X-axis. 

B.  Form  of  the  proposed  climate  index 

Because  the  climate  index  proposed  in  this  paper  is  an  indicator 
of  climatological  cloudiness,  it  is  instructive  to  review  the  status 
and  nature  of  current  global  cloud  climatologies.  As  mentioned  earlier 
good  cloud  climatologies  are  especially  needed  for  climate  modeling  and 
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climate  monitoring.  The  VIRES  climate  index  would  be  in  direct  compe¬ 
tition  with  standard  cloud  climatologies  to  fill  certain  of  these  re¬ 
quirements.  Therefore,  a  short  review  of  the  standard  approach  follows. 
A  more  detailed  synopsis  of  the  currently  available  cloud  climatologies 
is  given  in  two  reports,  one  by  Suomi  et  al.  (1977)  and  another  by 
Smith  (1978). 

The  most  widely  used  surface  based  observational  data  set  seems 
to  be  that  of  Telegados  and  London  (1954) .  Three  dimensional  distribu¬ 
tion  of  clouds  are  given  for  the  Northern  Hemisphere.  However,  the 
spatial  resolution  of  cloudiness  is  poor  and  as  a  consequence  even  the 
mean  zonal  cloudiness  values  given  are  of  questionable  accuracy. 

Climatologies  of  clouds  from  a  combination  of  ground  and  satellite 
sources  are  given  by  Sherr  et  al.  (1968)  and  by  Fye  (1978).  The  latter 
reference  describes  the  U.S.  Air  Force's  three  dimensional  nephanalysis 
model.  Both  sets  of  climatologies  suffer  from  the  lack  of  spatial  and 
temporal  homogeneity  in  the  quality  of  the  compiled  data.  Nevertheless 
both  climatologies  represent  global  coverage  and  have  three  dimensional 
cloud  fields.  Furthermore,  the  Air  Force's  data  base  is  continually 
being  added  to. 

The  longest  uniform  time-series  of  cloud  satellite  data  are  in  the 
form  of  visible  wavelength  brightness  values.  From  brightness  values 
total  cloud  amount  is  inferred.  Such  climatologies  are  given  by  Sadler 
(1969),  Sadler  et  al.  (1976),  Miller  and  Feddes  (1971),  and  Environ¬ 
mental  Satellite  Imagery  (1975-1976).  Steiner  (1978)  summarized  the 
Sadler  et  al.  data,  and  Avaste  et  al.  (1979)  attempted  to  compare  and 
combine  all  the  brightness  derived  cloud  climatologies  given  above. 


Avaste  found  that  the  main  limitation  of  this  data  is  its 


83 


non- comparability  and  lack  o£  accuracy.  There  is  also  no  cloud  height 
information.  This  last  drawback  almost  eliminates  climate  and  cloud 
modeling  applications. 

Inferences  of  clouds  from  satellite  data  using  other  methods  (see 
Chapter  II)  have  been  limited  to  brief  periods  and  usually  less  than 
global  coverage.  Although  most  methods  strive  for  three  dimensional 
information,  they  suffer  from  a  lack  of  accuracy  mainly  due  to  errors 
in  retrieval  technique  assumptions.  As  mentioned  earlier,  all  the 
approaches  use  standard  cloud  classification.  Some  reduce  the  classi¬ 
fication  to  3  heights  (low,  middle  and  high)  and  5  fractional  cate¬ 
gories,  while  others  use  more  elaborate  classifications  in  an  attempt 
to  be  more  descriptive.  Against  this  background  the  proposed  VIRES 
climate  index  will  be  discussed. 

Figures  17-23  show  the  final  composite  form  of  the  infrared  broad¬ 
band  weighting  functions  for  16  different  distributions  of  cloudiness 
for  the  typical  tropical  atmosphere.  Because  it  is  too  awkward  to  use 
the  whole  curve  for  purposes  of  an  index,  a  simpler  approach  is  needed. 
Of  course  the  two  curve  shape  parameters  define  the  curve  for  a  single 
case  situation,  but  they  will  not  define  the  curve  for  the  composite 
case.  Besides  being  simple,  the  index  should  be  unique  or  nearly  so 
for  a  given  cloud  distribution. 

After  much  searching  an  index  containing  four  parameters  was 
found  to  be  suitable.  The  first  parameter  of  the  index  represents  the 
level  of  the  atmosphere  above  which  25%  of  the  energy  emitted  to  space 
originates .  Likewise  the  next  two  parameters  represent  the  50%  and 
75%  levels.  The  fourth  number  represents  the  percentage  energy  escap¬ 
ing  to  space  from  the  surface.  The  approach  described  above  can  be 
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illustrated  using  the  following  equation: 


/  "*  THE?  “  lni> 

pi _ 


i. 


where  L  +  is  the  broadband  Infrared  radiance  measured  at  satellite 
sat 

level,  is  one  of  the  first  three  VIRES  index  numbers,  and  i  is  .25, 
.50  and  .75.  Of  course  other  i  values  may  be  chosen  to  define  another 
index. 

Table  14  gives  the  calculated  index  values  for  25  atmospheric 
VIRES  using  a  tropical  and  mid- latitude  profile  atmosphere.  The  atmos¬ 
pheric  levels  specified  in  the  index  are  in  mb.  For  reference  purposes 
Table  15  gives  the  conversion  from  mb  to  meters.  The  RTE  model  used 
produced  index  values  to  the  nearest  10  nd>.  Table  14  shows  no  surface 
emission  from  the  surface  in  the  tropical  case  due  to  water  vapor  dia¬ 
meter  absorption  and  approximations  used  in  the  RTE  model.  However, 
ideally  a  small  amount  of  surface  emission  would  be  expected. 

Index  values  given  in  Table  14  show  that  weighting  curves  that  are 
nearly  identical  give  identical  indexes  (i.e.  FLATS  and  FLATR,  BIMOS 
and  BIMOR,  1PK3S  and  1PK3R).  Specified  cloud  distributions  in  the  two 
different  atmospheres  give  the  same  index  at  high  levels  (i.e.  1PK3) 
but  a  somewhat  different  index  at  lower  atmospheric  levels  (i.e.  1PK8) 
due  primarily  to  the  fact  that  low  level  water  vapor  contributes  to  the 
atmospheric  VIRES.  Notice  that  low  clouds  (1PK7,  1PK8,  1PK9)  have 
nearly  the  same  index  in  the  tropical  case  where  low  level  water  vapor 


Pd  >M  fn 


Effective  Cloud 
Top  Distribution 


TROPCIAL  ATMOSPHERE  HID-LAT.  SUMMER  ATM. 
25Z  50Z  75Z  Sfc  Z  25Z  50Z  75Z  Sfc  Z 


FLATS 

250 

430 

650 

0 

250 

450 

690 

1 

FLATR 

250 

430 

650 

0 

250 

450 

690 

5 

FLATR* 

250 

430 

650 

0 

250 

450 

690 

6 

FLAT1 

250 

430 

650 

0 

250 

450 

690 

1 

BIMOS 

270 

410 

670 

0 

270 

430 

750 

0 

BIMOR 

270 

410 

670 

0 

270 

430 

750 

6 

BIM01 

230 

370 

690 

0 

230 

370 

730 

1 

BIM02 

290 

390 

630 

0 

290 

390 

630 

0 

BIM03 

270 

390 

550 

0 

270 

390 

550 

0 

1PK3S 

250 

350 

510 

0 

250 

350 

510 

0 

LPK3R 

250 

350 

510 

0 

250 

350 

510 

4 

1PK2 

190 

270 

390 

0 

190 

270 

390 

3 

1PK3 

230 

330 

430 

0 

230 

330 

430 

0  + 

1PK4 

270 

430 

510 

0 

270 

430 

510 

0 

1PK5 

290 

490 

590 

0 

310 

510 

590 

0 

1PK6 

310 

550 

670 

0 

310 

570 

690 

0 

1PK7 

310 

570 

730 

0 

330 

610 

770 

0  + 

1PK8 

310 

590 

750 

0 

330 

630 

850 

0  + 

1PK9 

310 

590 

750 

0 

330 

650 

910 

0  + 

CLEAR 

350 

650 

790 

0 

370 

790 

- 

37 

FLAT1 

290 

530 

730 

0 

310 

570 

950 

21 

BIM01 

270 

530 

730 

0 

270 

590 

930 

21 

1PK3 

290 

430 

710 

0 

290 

450 

930 

12  + 

1PK5 

310 

530 

710 

0 

330 

550 

890 

20 

1PK7 

330 

610 

750 

0 

330 

670 

890 

19 

+  Referred  to  in  the  text 


Table  14.  Proposed  VIRES  climate  index 


86 


Table  15.  Pressure  to  Height  Conversion. 


Pressure 

njb 

Mid.  Lat. 
Atm.  (m) 

Tropical 
Atm.  (m) 

110 

16280 

16330 

130 

15150 

15200 

150 

14150 

14230 

170 

13330 

13430 

190 

12600 

12710 

210 

11930 

12050 

230 

11340 

11450 

250 

10780 

10880 

270 

10260 

10360 

290 

9760 

9860 

310 

9290 

9390 

330 

8850 

8940 

350 

8430 

8510 

370 

8030 

8120 

390 

7640 

7720 

410 

7270 

7351 

430 

6920 

6990 

450 

6570 

6650 

470 

6240 

6320 

490 

5930 

6000 

510 

5620 

5690 

530 

5320 

5390 

550 

5030 

5090 

570 

4750 

4810 

590 

4480 

4540 

610 

4210 

4270 

630 

3960 

4010 

650 

3710 

3750 

670 

3460 

3510 

690 

3220 

3260 

710 

3000 

3030 

730 

2760 

2800 

750 

2540 

2570 

770 

2320 

2350 

790 

2110 

2130 

810 

1900 

1920 

830 

1690 

1720 

850 

1490 

1520 

870 

1300 

1320 

890 

1110 

1120 

910 

917 

933 

930 

732 

745 

950 

550 

561 

970 

372 

379 

990 

197 

201 

1007 

55 

56 
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emits  like  a  cloud.  This  is  not  true  for  the  mid-latitude  case  where 
low  level  water  vapor  does  not  contribute  much  to  the  atmospheric 
VIRES.  Also  notice  that  similar  distributions  of  clouds  but  with  a 
different  percentage  of  clear  sky  have  different  indexes  (l.e.  clear 
1PK3  vs.  1PK3) . 

In  general  the  VIRES  climate  index  outlined  above  may  be  consider¬ 
ed  an  appropriate  and  accurate  atmospheric  descriptor.  The  parameters 
in  the  index  indicate  the  vertical  source  of  emitted  earth-atmospheric 
energy.  The  index  is  especially  compatible  with  the  original  satel¬ 
lite  data  resulting  in  minimal  anticipated  errors.  The  analysis  shows 
the  index  has  the  desired  uniqueness  characteristic  while  at  the  same 
time  remaining  straightforwardly  simple. 

Of  course  one  could  consider  it  a  drawback  that  the  index  is  not 
compatible  with  standard  measures  of  cloudiness.  However,  the  index 
is  meant  to  be  closely  related  to  the  radiative  measurements.  Clouds 
are  thought  of  as  simply  one  constituent  (albeit  the  most  Important) 
of  the  atmospheric  VIRES.  Furthermore  real  physical  clouds  are  never 
considered  at  all.  Only  the  equivalent  radiative  effects  of  specified 
clouds  are  taken  into  account.  This  approach  avoids  many  difficulties 
while  at  the  same  time  providing  information  that  may  be  used  for 
climate  monitoring  and  for  calibration  and  parameterization  of  radia¬ 
tive  calculations  in  climate  models  (the  topic  of  Chapter  VII) .  The 
next  section  discusses  possible  specifications  for  index  use. 

C.  Climate  index  specifications 

The  climate  index  described  in  this  paper  will  be  used  similar  to 
a  cloud  climatology.  It  should,  therefore,  have  some  of  the  same 
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characteristics.  Smith  (1978)  suggests  that  for  purposes  of  param¬ 
eterization  of  cloudiness  and  radiation  in  climate  models,  cloud 
climatologies  should  have  global  coverage  with  a  space  resolution  of 
250  km  (s  2*s°  Lat.  by  2%#  Long.)  and  a  time  resolution  of  4  hours.  Re 
search  by  Avaste  et  al.  (1979)  suggests  that  for  climate  monitoring 
purposes  the  space  resolution  be  500  km  (~  5s  Lat.  by  5°  Long.)  and 
time  resolution  be  monthly.  A  weekly  specified  time  resolution  has 
been  suggested  by  the  U.S.  Committee  for  GARP  (1975).  Before  specify¬ 
ing  time  and  space  resolution  for  the  proposed  climate  index  let  us 
consider  satellite  data  collection  limitations. 

Present  operational  and  experimental  sun-synchronous  weather 
satellites  carry  infrared  radiometers  that  have  nadir  resolutions  of 
at  least  30  km.  Similar  radiometers  for  geosynchronous  satellites 
will  have  comparable  resolution.  Therefore,  compiling  a  VIRES  index 
climatology  with  useful  space  resolution  seems  feasible.  However,  the 
global  coverage  requirement  Is  harder  to  meet  when  combined  with  the 
time  resolution  demand. 

Satellites  such  as  NIMBUS  and  TIROS  are  in  sun-synchronous  orbits 
and  only  make  observations  over  any  single  earth  location  twice  a  day 
at  the  same  times  everyday  although  the  whole  globe  is  covered.  Be¬ 
cause  of  diurnal  cloudiness  variations,  twice  a  day  observations  give 
biased  results  (Harrison  et  al.  1978,  1980).  On  the  other  hand,  geo¬ 
synchronous  satellites  maUe  observations  at  all  hours  but  in  general 
can  effectively  view  only  between  60°N  and  60°S  latitude  and  about  120 
of  longitude.  More  complete  coverage  can  be  achieved  by  observing 
nearly  all  the  low  and  middle  latitudes  at  each  hour  by  adding  a  satel 
lite  with  a  mid-inclined  orbit  which  precesses  (Harrison  et  al.  1978). 


A  good  discussion  concerning  optimizing  satellite  observations  is  given 
by  Campbell  and  Vonder  Haar  (1978).  Nevertheless,  it  is  clear  that  no 
one  type  of  satellite  offers  both  the  ideal  spatial  coverage  and  ideal 
temporal  resolution. 

Remembering  this  limitation,  one  should  consider  the  importance  of 
cloud  classification.  Sherr  et  al.  (1968)  proposed  29  different  cloud 
climatological  regions  globally.  These  regions  were  based  primarily 
on  seasonal  distributions  of  mean  monthly  cloud  cover.  Also  considered 
were  annual  cloud  distributions,  area  precipitation  distributions,  and 
different  climate  classification  schemes.  Most  regions  were  repeated 
two  or  more  times  throughout  the  world.  Some  of  Sherr' s  typical  region 
descriptions  are;  tropical  cloudy,  desert  marine,  mid-latitude  clear 
summer,  high  latitude  clear  winter,  mediterranean  and  polar.  A  simi¬ 
lar  cloud  classification  scheme  is  given  by  Winston  (1969). 

By  considering  cloudiness  regions  instead  of  strictly  global 
measurements  one  might  for  instance  actually  have  a  more  sensitive 
measure  for  the  purpose  of  climate  monitoring.  Applications  of  the 
climate  modeling  type  might  also  be  stronger  on  a  regional  basis. 
Consequently,  global  coverage  seems  to  be  of  secondary  importance. 

Thus,  one  would  expect  a  climate  index  derived  from  geosynchronous 
satellite^ data  to  be  the  basis  of  a  good  climatology  except  where  total 
global  coverage  is  required.  Sun-synchronous  satellites  could  provide 
additional  information  if  the  time  bias  problem  is  avoided. 

From  this  discussion  one  may  conclude  that  the  best  approach  is 
to  take  the  geostationary  satellite  information  available  and  average 
it  to  a  250  km  space  resolution  and  4  hour  time  resolution.  From 
these  values  longer  time  and  larger  space  scale  values  can  be  obtained 
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with  emphasis  on  weekly  regional  averages.  Of  course  there  should  be 
an  attempt  to  use  sun-synchronous  satellite  data  to  complete  the  global 
picture . 

Related  to  the  time  resolution  problem  is  the  problem  of  instru¬ 
ment  continunity  through  time.  In  other  words,  how  do  we  keep  the 
satellite  instrument  absolutely  calibrated?  Without  the  absolute  cal¬ 
ibration  it  becomes  necessary  to  intercompare  instruments  from  dif¬ 
ferent  satellites  for  the  purpose  of  obtaining  a  homogeneous  data  set. 
It  is  also  necessary  to  detect  and  correct  changes  in  a  specific  in¬ 
strument  over  time.  Of  course  intercomparisons  would  be  useful  as  a 
second  check  even  with  absolute  calibration.  The  proposed  Space 
Shuttle  may  offer  the  means  of  intercomparing  instruments,  although 
for  now  such  a  process  is  impossible  (COSPAR,  1978b).  Also,  a  common 
ground  base  laboratory  test  facility  for  all  satellite  instruments 
would  provide  a  chance  to  uncover  systematic  differences. 

Calibrating  radiation  instruments  in  any  absolute  sense  is  very 
difficult  if  not  impossible  (COSPAR,  1978b)  especially  for  SW  sensors. 
Calibration  targets  such  as  black  bodies  are  sometimes  used  success¬ 
fully  for  IR  instruments .  One  may  also  use  a  transfer  technique;  ob¬ 
serve  the  same  target  at  the  same  time  with  two  different  instruments 
whose  output  can  be  compared.  If  one  of  the  instruments  has  a  trusted 
calibration  then  the  other  one  may  be  calibrated.  Another  approach  is 
to  convert  radiant  energy  into  heat  energy  using  a  known,  stable  pro¬ 
cess.  The  advantage  here  is  that  the  efficiency  of  this  conversion  as 
a  function  of  temperature  can  be  determined  independently  and  therefore 
'known'  calibration  sources  are  not  needed.  Lastly,  it  can  be  said 
that  the  nature  of  radiation  (e.g.  its  amplitude,  wavelength. 
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Interference  and  diffraction  effects,  and  polarization)  makes  calibra¬ 
tion  of  radiometers  a  complex  operation  with  many  uncertainties,  es¬ 
pecially  for  SW  instruments.  Fortunately,  the  more  easily  calibrated 
IR  sensitive  radiometers  are  used  to  derive  the  VIRHS  index. 
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VII.  POSSIBLE  APPLICATIONS  OF  A  SATELLITE  DERIVED  CLOUDINESS  INDEX  j 

1 

i 

t 

Most  of  the  cloud  retrieval  methods  mentioned  in  Chapter  II  have  j 

been  developed  with  the  idea  that  deducing  global  cloud  climatology  in¬ 
formation  from  satellite  data  is  desirable.  Most  have  tried  to  conform 
with  classical  ground  observational  definitions  of  clouds.  The  present 
research  suggests  that  a  more  realistic  and  progressive  point  of  view 
is  to  consider  clouds  in  light  of  their  radiative  properties.  After 
all,  counting  and  typing  individual  clouds  on  a  global  climatological 
basis  is  probably  an  impossible  task,  and  the  results  of  such  a  census 
may  not  provide  the  best  information  for  the  application.  Therefore,  a 
climatic  index  that  represents  cloudiness  by  means  of  the  infrared 
broadband  weighting  function  is  proposed.  How  might  such  an  index  be 
used?  Perhaps  the  two  most  important  applications  for  such  a  climato¬ 
logical  index  are  in  the  areas  of  climate  modelling  and  climate  moni¬ 
toring  (Smith,  1978). 

A  quote  from  GARP  Publication  16  (Stockholm,  1974)  relates  the 
importance  between  climatic  processes  of  clouds  and  their  radiative 
effects:  "...  proper  treatment  of  radiative  effects  of  clouds  is  the 
single  most  important  factor  in  the  overall  parameterization  of  radia¬ 
tion  in  a  climate  model".  However,  proper  treatment  of  clouds  is 
difficult  even  for  the  most  advanced  climate  models.  For  example,  most 
clouds  are  sub-grid  scale  for  the  typical  general  circulation  model 
(GCM).  Thus,  clouds  will  undoubtedly  be  parameterized  in  future  models 
much  as  they  are  now.  Furthermore,  these  cloud  parameterizations  are 
one  of  the  weakest  aspects  of  the  present  GCM's.  Climate  models  now 
lack  the  ability  to  meaningfully  calculate  radiative  forcing  due  to  the 
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failure  of  cloud- dynamic  parameterizations  (Webster,  1978).  Neverthe¬ 
less,  as  newer  and  better  cloud  models  are  developed  they  must  be  veri¬ 
fied  and  tuned  against  actual  measurements  of  global  cloudiness  (Cox, 
1978). 

Because  climate  models  necessarily  have  many  assumptions,  approxi¬ 
mations  and  parameterizations,  verification  of  their  results  is  quite 
important.  Furthermore,  computation  of  one  parameter  may  be  successful 
because  the  model  contains  compensating  errors  in,  for  example,  the 
treatment  of  physical  processes.  One  can  have  more  confidence  in  the 
validity  of  ones  model  by  checking  calculated  against  measured  param¬ 
eters  that  are  both  a  direct  measure  of  a  physical  process  and  avail¬ 
able  in  time  and  space  scale  detail.  The  climatic  cloudiness  index 
outlined  in  this  paper  is  a  particularly  useful  parameter  because  it  is 
a  radiative  index  that  can  be  compiled  for  weekly  or  seasonal  world¬ 
wide  vertical  profile  values.  In  other  words,  the  4-dimensional  char¬ 
acteristic  of  this  observable  variable  would  make  it  a  particularly 
good  verification  tool  for  the  parameterization  of  clouds  in  the  GCM 
radiative  calculations . 

The  values  of  the  proposed  index  may  also  be  useful  in  GCM  cloud 
parameterization  studies.  Instead  of  dealing  directly  with  clouds  in 
a  model,  one  might  go  directly  from  dynamic-thermodynamic  considera¬ 
tions  to  radiative  effects  (COSPAR,  1978a) .  Another  approach  is  to 
relate  model  derived  cloudiness  directly  to  radiative  divergence  pro¬ 
files  without  any  model  cloud  radiative  calculations  (Cox  and 
Vonder  Haar,  1973).  Such  "direct"  parameterization  studies  will  need 
to  use  observational  data  like  the  cloudiness  and  radiative  informa¬ 
tion  contained  in  the  index. 
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The  second  area  of  potential  index  application  is  in  climate  moni¬ 
toring.  The  earth's  climate  is  a  function  of  the  earth-atmosphere  sys¬ 
tem  reacting  to  the  equator- to-pole  gradient  of  net  energy;  this  net 
energy  budget  may  be  quantified  and  divided  into  components  by  the  net 
radiation  budget  studies  using  satellite  data  (Vonder  Haar  and  Oort, 
1973;  Campbell  and  Vonder  Haar,  1980a, b).  Clouds  are  the  principal 
modulators  of  the  radiation  for  they  strongly  Influence  the  earth's 
solar  albedo  and  infra-red  absorption.  To  extend  and  Improve  the  cli¬ 
mate  monitoring  aspect  of  these  radiation  budget  studies  there  is  a 
need  to  determine  cloudiness  independent  of  the  satellite  radiation 
budget  measurements  (Vonder  Haar,  1979).  The  index  described  in  this 
paper  would  be  useful  for  this  purpose. 

Because  we  feel  that  cloudiness  is  an  important  climate  variable, 
and  since  the  VIRES  index  is  strongly  a  function  of  cloudiness,  it 
follows  that  the  index  itself  represents  a  physical  characteristic  of 
the  climate  system.  The  temporal  and  spatial  variations  of  this  index 
can  be  detected  using  the  statistical  method  described  earlier.  This 
could  be  done  on  a  global  or  regional  scale  in  an  attempt  to  detect 
climatic  trends.  Regional  changes  which  can  be  concealed  in  global 
averages  may  provide  a  more  sensitive  measure  of  climatic  change. 
Moreover,  climatic  anomalies  in  one  area  are  often  correlated  with 
variations  in  another  area.  Occassionally  there  is  also  a  time  lag. 

In  such  cases  climatic  forecasts  can  be  made.  These  teleconnections 
might  be  identified  and  observed  easier  by  processing  satellite  data 
into  meteorological  Indexes  like  the  one  described  in  this  paper. 

The  index  described  in  this  paper  has  both  information  about 
cloudiness  and  outgoing  longwave  radiation.  In  this  respect  it  would 
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be  especially  valuable  to  the  National  Climatic  Research  Program 
(NCRP) .  The  NCRP  Committee  proposes  satellite  observations  in  the  form 
of  a  Climatic  Index  Monitoring  Program  in  addition  to  conventional 
meteorological  observations  (U.S.  Committee  for  GASP,  1975).  This 
VIRES  index  would  represent  two  indexes  in  one.  Furthermore,  it  can  be 
computed  for  the  global  coverage  and  weekly  frequency  required. 

There  are  other  possible  uses  for  a  radiatlvely  defined  cloudiness 
index.  When  combined  with  radiation  budget  studies  a  better  assessment 
of  atmospheric  energetics  may  be  possible  (COSPAR  WG  6,  1978a).  Cli¬ 
matic  change  may  be  viewed  as  an  adjustment  among  compensating  feedback 
processes.  Using  this  index  one  may  be  able  to  Identify  and  quantify 
coupling  between  cloudiness  and  other  atmospheric  variables.  Perhaps 
by  using  the  index  more  specific  relationships  between  the  cloudiness 
and  the  earth-atmosphere  energy  balance  (Hartman  and  Short,  1980; 

Herman  et  al.  1980;  Ohring  and  Clapp,  1980)  can  be  found.  Conceivably 
this  index  could  help  determine  whether  or  not  a  change  in  cloudiness 
necessarily  results  in  a  change  in  the  climate  (Cess,  1976;  Ellis, 

1978;  Coakley,  1979;  Van  Den  Dool,  1980).  In  this  regard  the  index 
could  shed  some  much  needed  light  on  the  role  of  cloud  vertical  struc¬ 
ture  as  it  applies  to  the  cooling  to  space  portion  of  the  earth- 
atmosphere  radiation  budget.  The  VIRES  index  would  also  be  useful  in 
understanding  more  clearly  the  consequence  of  tropospheric  water  vapor 
emission  as  it  relates  to  the  earth-atmosphere  radiative  loss. 
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VIII.  CONCLUSIONS 

This  paper  introduces  and  describes  a  climate  index  called  the 
VIRES  index.  VIRES  is  the  acronym  for  Vertical  Infrared  Radiative 
Emitting  Structure  and  quite  by  chance  it  is  also  the  plural  form  of 
the  Latin  word  vis,  meaning  forces  or  powers.  Thus,  the  acronym  seems 
especially  appropriate  since  the  atmosphere's  VIRES  is  one  of  the  major 
forcing  factors  behind  the  earth's  climate.  The  VIRES  index  is  related 
to  the  earth's  climate  through  the  earth's  radiation  budget  and  there¬ 
fore  can  be  considered  a  climate  index.  The  logic  behind  this  claim 
may  be  stated  as  follows.  The  climate  system  is  determined  by  the 
energy  input  to  the  system  and  the  distribution,  transformation  and 
storage  of  energy  in  various  forms  within  the  system.  These  processes 
are  mirrored  in  the  components  of  the  earth's  radiation  budget,  one  of 
which  is  the  outgoing  emitted  thermal  radiation  (COSPAR  Report  to  ICSU 
and  JOC,  1978b).  This  cooling  to  space  is  described  by  the  VIRES  which 
is  primarily  a  function  of  the  cloud  distribution. 

This  index  is  an  attempt  to  optimize  the  use  of  satellite  data 
for  climate  purposes  by  directly  utilizing  the  radiative  aspects  of 
the  atmosphere  while  avoiding  some  of  the  difficulties  of  inferring 
standard  meteorological  variables  from  satellite  radiances.  The  VIRES 
index  is  based  on  broadband  infrared  weighting  curves  retrieved  from 
operationally  measured  spectral,  earth-emitted  radiation  in  the  CO2 
absorption  band.  These  curves  describe  the  vertical  structure  of  in¬ 
frared  radiative  emission  and  are  a  function  of  the  cloud,  temperature 
and  moisture  distributions.  The  most  important  findings  of  this  re¬ 
search  are  summarized  below. 
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A.  Cloud  -  VIRES  relationship 

This  paper  demonstrates  under  both  tropical  and  mid-latitude 
atmospheric  conditions  the  predominate  influence  of  the  three  dimen¬ 
sional  cloud  distribution  on  the  atmosphere's  Vertical  Infrared  Radia¬ 
tive  Emitting  Structure  (VIRES).  Therefore,  the  VIRES  index  can  be 
considered  an  indicator  of  cloudiness  with  the  caution  that  low  cloud 
or  high  water  vapor  content  can  both  result  in  the  same  VIRES.  Also 
illustrated  is  the  fact  that  the  atmospheric  VIRES  is  responsive  to 
the  three  dimensional  distributions  of  water  vapor  and  temperature. 

B.  VIRES  retrieval  technique 

VIRES  curves  may  be  Inferred  directly  from  satellite  measured 
radiances.  A  retrieval  technique  utilizing  CO^  band  spectral  radiances 
is  described  which  determines  two  weighting  curve  shape  parameters 
(pwj  and  a) .  The  empirical  RTE  method  was  found  to  be  both  computa¬ 
tionally  fast  and  accurate.  The  two  parameters  inferred  using  the 
empirical  RTE  method  define  a  single  scene  VIRES.  The  single  scene 
retrievals  are  composited  to  give  a  temporal  and  spatial  average  VIRES 
curve.  A  sensitivity  study  and  error  analysis  using  simulated  satel¬ 
lite  data  that  included  the  effect  of  sampling  inadequacies  quantified 
the  abilities  of  the  empirical  RTE  retrieval  approach. 

C.  VIRES  uniqueness 

The  VIRES  curves  calculated  for  diverse  cloud  (climate)  regimes 
have  been  shown  to  be  statistically  different  in  all  cases  with  the 
following  two  exceptions.  In  the  moist  atmosphere  case  the  high  con¬ 
centrations  of  low  level  water  vapor  radiatively  emit  to  space  almost 
identically  like  the  case  with  low  level  cloud.  The  second  exception 
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occurs  for  comparisons  between  two  similar  cloud  shape  distributions. 
The  simple  non-parametrlc  statistical  test  compares  the  points  from 
two  VIRES  curves  every  20  mb.  As  might  be  expected,  the  VIRES  index 
values  show  (qualitative)  uniqueness.  Thus,  climate  change  can  be 
effectively  monitored  using  the  VIRES  index  which  conveniently  and 
quantitatively  expresses  the  vertical  structure  of  terrestrial  emission 
to  space. 


D.  VIRES  index 

Although  the  VIRES  curves  contain  the  maximum  vertical  informa¬ 
tion,  they  are  cumbersome  to  work  with.  Therefore,  a  useful  descriptor 
of  the  VIRES  curve  called  the  VIRES  index  was  developed.  As  stated 
above,  the  index  is  able  to  protray  differences  in  the  atmosphere’s 
vertical  emitting  structure.  The  first  three  number  parameters  in  the 
index  represent  pressure  levels  in  the  atmosphere  above  which  a  speci¬ 
fied  fraction  (.25,  .50,  .75)  of  the  infrared  energy  lost  to  space 
originates.  The  last  parameter  represents  the  fraction  of  energy  lost 
to  space  which  originated  at  the  earth's  surface.  The  index  is  es¬ 
pecially  useful  since  it  represents  the  VIRES  curve  information  in  a 
shortened,  interpre table,  and  flexible  form.  The  four  parameter  in¬ 
dex  prposed  in  this  work  may  be  easily  modified  to  maximize  its 
usefulness  for  different  applications. 

E.  Suggested  applications  and  sampling  strategies 

A  number  of  specific  applications  have  been  proposed  and  in 
general  they  are  related  to  climate  monitoring  and  climate  modeling. 

For  example,  to  monitor  the  climate,  regional  VIRES  index  averages  can 
be  compiled.  These  values  would  probably  be  very  sensitive  to  climate 


variations.  It  is  recommended  that  this  basic  index,  calculated  from 
geostationary  satellite  radiance  measurements,  be  averaged  to  a  spatial 
resolution  of  250  km  (2  1/2®  latitude  by  2  1/2®  longitude)  with  the 
smallest  time  resolution  being  4  hours.  Of  course  sun  synchronous 
satellite  data  may  be  used  to  supplement  coverage.  From  these  data 
coarser  spatial  and  temporal  averages  can  be  calculated  for  other 
specific  applications  such  as  regional  or  global  climate  modeling.  The 
VIRES  index  can  be  used  for  climate  model  baselining  and  verification 
of  radiative  calculations.  Index  values  may  also  be  useful  to  research 
efforts  in  the  area  of  cloud  parameterlzations. 
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APPENDIX 

Retrieval  Error  Analysis 

Chapter  V  Is  an  Introduction  and  summary  of  the  analysis  presented 
below.  For  an  explanation  of  the  boundary  conditions  and  assumed  error 
values  refer  to  Chapter  V.  This  Appendix  contains  the  following  de¬ 
tailed  error  analysis.  Errors  due  to  system- instrument  errors  are 
investigated  first,  followed  by  errors  due  to  temperature  and  humidity 
profile  unknowns.  The  magnitude  of  errors  due  to  Ct^  profile  varia¬ 
tions  are  also  examined.  Errors  due  to  specified  cloud  properties  are 
reviewed,  followed  by  a  brief  analysis  of  combined  error  sources.  The 
last  source  of  error  examined  is  the  assumption  that  a  single  repre¬ 
sentative  radiating  surface  is  in  the  radiometer  field  of  view.  The 
following  sections  detail  the  error  findings. 

A.  Random  sensor  error  effects 

Figures  24-31  give  the  bias  and  RMS  retrieval  errors  due  to  random 
sensor  errors  for  both  a  mid-latitude  and  tropical  atmosphere.  In¬ 
strument  system  noise  is  assumed  Gaussian  with  a  mean  of  zero  and  a 

-2  -1 

standard  deviation  of  0.22  raW  m  sr  cm  for  the  15  ym  channels,  and 
-2  -1 

0.11  mW  m  sr  cm  for  the  window  channel.  These  values  represent  the 
state  of  the  art  precision  of  the  HIRS  instrument  on  TIROS-N  (Schwalb, 
1978).  For  calculation  purposes  the  set  of  spectral  values  represent¬ 
ing  a  particular  atmosphere  are  modified  with  random  noise  errors  that 
are  distributed  as  described  above.  These  modified  spectral  radiance 
values  are  then  used  in  the  empirical  RTE  curve  shape  parameter  re¬ 
trieval  algorithm  to  solve  for  the  calculated  pw^  and  a  values,  which 
are  then  compared  to  the  known  correct  values.  Statistics  are 
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Figure  24.  pw^  pressure  level  bias:  random  sensor  errors  for  mid¬ 
latitude  summer  model. 
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Figure  28.  a  bias:  Random  sensor  error  for  mid-latitude  summer  model 
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Figure  29.  a  RMS:  Random  sensor  error  for  mid-latitude  summer  model 
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calculated  using  60  sets  of  calculated  versus  given  values  for  each 
particular  atmospheric  condition.  As  seen  In  Figures  24-31,  one 
hundred  statistics  per  figure  are  calculated  using  combinations  of  10 
effective  cloud  top  pressures  and  10  effective  cloud  fraction  values. 

Understanding  the  bias  statistic  in  Figures  24,  26,  28  and  30  is 
simplified  if  one  realizes  that  there  is  a  correlation  between  the  re¬ 
trieved  values  of  pw^  and  a.  For  example,  a  single  spectral  radiance 
may  be  the  result  of  small  pw^  and  a  values  or  the  result  of  large  pw^ 
and  a  values.  For  the  single  spectral  radiance  there  is  a  spectrum  of 

p  ,  and  a  value  pairs  going  from  both  small  to  both  large  that  could 
wi 

define  a  weighting  curve  for  that  radiance  value.  Thus,  to  satisfy  a 
set  of  spectral  radiances  with  instrument  noise  the  empirical  routine 
will  usually  tend  to  err  in  the  same  direction  for  both  pwf  and  a.  For 
the  most  part.  Figures  24  and  28,  and  Figures  26  and  30  show  this  error 
characteristic . 

Furthermore,  the  bias  figures  show  errors  due  principally  to 
boundary  conditions.  The  physical  constraints  of  the  system  force  the 
instrument  noise  error  to  show  bias  errors  when  pw^  -*■  1000  mb.  In  this 
case  radiating  surfaces  below  1000  mb  are  not  allowed  causing  the  p^ 
bias  error  to  be  negative  (bias  *  calculated  -  given) .  At  this  boun- 
ary  a  bias  errors  are  also  generally  negative.  Figures  24  and  28  show 
the  same  effect  as  a  -*■  0.  By  necessity  the  bias  is  toward  larger  a 
values  and  as  a  consequence  larger  pw£  values  (positive  bias).  Figures 
26  and  30  show  this  effect  to  a  lesser  degree.  Overall  this  analysis 
shows  the  retrieval  is  highly  unbiased  due  to  instrument  noise  except 
as  a  -*■  0  and  p  ,  -*■  1000  mb.  As  stated  previously,  when  a  is  near  zero 
and  pw£  is  near  the  surface  the  broadband  weighting  curve  is  only 
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slightly  changed  from  the  clear  sky  case.  Thus,  these  error  effects 
are  minimized. 

Figures  25,  27,  29  and  31  give  the  RMS  error  due  to  system  noise 
of  pwj  and  a  for  a  mid-latitude  and  tropical  atmosphere.  The  more 
gradual  increase  in  RMS  errors  for  pw^  and  a  as  a  +  0  and  p  ^  -*■  1000 
mb  indicates  that  the  RMS  error  is  probably  related  closely  to  signal 
to  noise  ratio  which  gradually  decreases  as  a  +  0  and  p  ^  •+■  1000  mb. 

One  also  notices  that  RMS  values  are  slightly  larger  for  the  mid¬ 
latitude  case  (Figures  25  and  29).  This  is  probably  due  to  the  fact 
that  the  tropical  atmosphere  contains  much  more  water  vapor.  The  pre¬ 
sence  of  water  vapor  may  act  as  a  buffer  reducing  the  impact  radiance 
errors  have  on  the  retrieval  process.  In  conclusion,  it  appears  that 
random  noise  effects  produce  py^  and  a  retrieval  errors  that  are  in  the 
mean  quite  small.  Any  specific  case  retrieval  errors  will  be  of  the 
magnitude  indicated  by  the  RMS  figures. 

B.  Profile  error  effects 

Analysis  of  profile  error  effects  is  extremely  difficult  due  to 
the  computational  time  involved.  Instead  of  adding  random  noise  60 
times  to  sets  of  spectral  radiance  values  (one  for  mid- latitude  and 
another  for  tropical  atmospheres) ,  the  sets  of  spectral  radiance  values 
have  to  be  generated  using  profiles  of  temperature,  humidity  and  CO^ 
that  are  modified  with  errors.  To  run  the  spectral  RTE  for  each  wave- 
number  for  60  different  error  plagued  profiles  for  each  of  the  two 
basic  atmospheres  in  an  attempt  to  produce  figures  similar  to  Figures 
24-31  requires  a  large  amount  of  computer  time.  Instead  of  taking  this 
route,  a  less  descriptive  but  much  faster  approach  is  used  and 
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outlined  below.  It  is  the  same  technique  described  at  the  end  of 
Chapter  IV. 

Bias  ar.d  RMS  errors  are  calculated  from  50  values.  These  values 
come  from  calculations  made  at  the  5  effective  cloud  top  levels  within 
a  100's  range  (i.e.  200,  220,  240,  260  and  280  mb)  for  10  fractional 
effective  cloud  amounts  (i.e.  1.0,  0.9,  ...»  0.1).  The  results  given 
are  only  for  the  mid-latitude  atmosphere  although  they  are  also  re¬ 
presentative  of  the  tropical  atmosphere  case. 

For  example.  Figure  32  which  represents  this  approach  using  in¬ 
strument  noise,  may  be  comprared  to  Figures  24,  25,  28  and  29.  Notice 
that  if  the  lines  in  Figure  32  are  considered  trends  representing  a 
single  overall  statistic,  with  no  importance  placed  on  variations  be¬ 
tween  effective  cloud  top  categories,  they  compare  favorably  with  the 
trends  indicated  in  Figures  24,  25,  28  and  29.  Thus,  the  trends  in¬ 
dicated  in  Figures  32-37,  especially  those  trends  neglecting  the  two 
end  groupings  are  taken  as  representative  of  a  mean  value  of  bias  or 
RMS.  Remenfcer  that  either  end  effective  cloud  category  shows  effects 
due  to  boundary  conditions.  Despite  obvious  limitations,  these  trend 
values  are  regarded  as  sufficiently  descriptive  for  this  analysis. 

Figure  33  shows  the  bias  and  RMS  of  a  and  pw^  associated  with  a 
2°C  RMS  and  50%  water  vapor  profile  uncertainty  specified  randomly  and 
independently  for  each  20  mb  layer.  The  distribution  of  errors  is 
Gaussian  with  zero  mean.  Water  mixing  ratio  errors  are  given  as  a 
percent  of  the  correct  mixing  ratio  at  any  given  level.  For  this 
particular  analysis  instrument  noise  is  set  to  zero. 

For  each  plot  in  Figure  33  there  are  two  curves.  The  dashed  line 
represents  the  error  when  the  specified  climatological  temperature  and 


(mb)  Pwf  RMS  (mb) 
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Effective  Cloud  Top  Cotegories 


Figure  32.  Random  sensor  error  effect  on  shape  parameter  retrieval. 


Figure  33.  +2°  RMS  temperature  and  50%  humidity  error  effects  on 

shape  parameter  retrieval. 


v  Bias  +  N  RMS  Bio*  (mb 

80S  8  o _ o  8  S  o _  8  8  5  00 


Effective  Cloud  Top  Categories 


36.  +  50X  LWC  error  effect  on  shape  parameter  retrieval 


humidity  profile  Is  used.  The  dotted  line  represents  error  analysis 
when  a  clear  column  radiance  correction  (CCRC)  is  applied.  As  can  be 
seen  in  Figure  33,  the  application  of  CCRC,  which  simply  adds  one  more 
piece  of  information  (^Vca)»  greatly  reduces  the  RMS  errors  and  seems 
to  improve  the  bias  error.  The  inclusion  of  the  measured  LVcg  informa¬ 
tion  is  reasonable  since  in  many  cases  this  information  will  be  avail¬ 
able.  Even  in  overcast  situations,  using  the  closest  L  ,  measurement 

vcs 

should  be  an  improvement  over  assuming  the  Lvcg  inferred  from  the 
climatological  profile  is  correct. 

The  CCRC  is  applied  in  the  following  way.  Let  L  be  the  clear 

sky  spectral  radiance  from  a  climatological  profile.  Let  Lcgm  be  the 

\ 

measured  clear  sky  spectral  radiance .  Each  spectral  radiance  value 

used  in  the  retrieval  procedure  is  then  multiplied  by  a  factor  equal  to 

CL  f  L  )  . 
esc  csm  V 

Figure  34  gives  the  RMS  bias  errors  due  solely  to  a  +  1%  error  in 
the  CO2  profile.  The  error  is  not  applied  randomly  but  through  the 
whole  profile  since  C02  is  well  mixed  in  the  troposphere.  The  magni- 
tude  represents  of  the  seasonal  maximum  change  in  C02  observed. 

Again  the  CCRC  technique  eliminates  the  effect  of  the  error  almost  en¬ 
tirely.  The  main  effect  of  the  C02  error  is  to  change  the  pwf  value 
from  true  while  a  is  calculated  correctly.  When  there  is  more  C02  than 
assumed  the  radiation  surface  is  retrieved  higher  in  the  atmosphere 
that  it  really  is.  Notice  a  change  of  scale  in  the  RMS  plots.  Even 
without  the  CCRC  applied  to  the  C02  error  the  temperature  and  humidity 
as  well  as  instrument  errors  are  far  more  important  in  terms  of  RMS 


error. 


C.  Effective  cloud  radiative  properties  effect 

As  i- tated  previously,  because  the  empirical  RTE  retrieval  method 
solves  for  an  a  value,  theoretically  the  emittance  for  the  radiating 
surface  being  retrieved  may  be  arbitrarily  specified.  The  accuracy  of 
this  statement  is  tested  in  this  section.  Errors  due  to  both  over  and 
under  specifying  emittance  will  be  examined  by  neglecting  all  other 
error  sources . 

Figure  35  shows  the  effect  of  assuming  effective  cloud  depths  of 
100  nb  when  they  are  actually  20  mb  thick.  In  other  words  the  emit¬ 
tance  is  assumed  to  be  much  larger  than  it  really  is.  By  comparing 
Figure  35  to  Figures  13  and  14  one  can  see  the  BMS  and  bias  errors  for 
Pw£  are  essentially  alike.  The  same  is  true  for  a  errors  which  are  not 
plotted.  Instead  the  errors  for  N  are  plotted  to  illustrate  that  if 
emittance  is  assumed  too  large  the  N  value  is  simply  decreased  to  give 
the  correct  a  value  (a  -  eN) .  Notice  that  below  600  mb,  20  mb  thick 
effective  clouds  behave  much  like  100  mb  thick  ones. 

Figure  36  gives  an  analysis  of  the  situation  where  the  emittance 
is  actually  greater  than  specified.  This  error  occurs  only  at  higher 
levels  since  emittance  below  300  mb  is  specified  at  its  maximum  value. 
Again  N  statistics  are  plotted  instead  of  a  statistics.  Since  emit¬ 
tance  is  assumed  too  small  at  the  upper  tropospheric  levels,  the  N 


value  calculations  are  biased  positive  to  give  the  correct  a  value. 

In  this  case  the  curves  show  small  variations  from  those  plotted  in 
Figures  13  and  14. 

Figure  37  shows  the  results  of  combining  errors.  The  errors 
used  are  a  random  2°C  RMS  temperature  profile  error  with  50%  variation 
in  humidity  as  explained  in  Section  B.  Also  included  is  random 
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Instrument  noise  as  explained  In  Section  A.  Furthermore,  true  cloud 
depth  Is  allowed  to  vary  randomly  from  20  to  180  mb  and  LWC  is  allowed 
to  vary  randomly  f rom  +  50%  of  the  specified  value. 

The  consistency  of  the  RMS  error  from  one  effective  cloud  top 
category  to  another  suggests  that  the  limited  sample  analyzed  gives  a 
typical  value  of  RMS.  However,  confidence  in  the  bias  errors  cannot  be 
as  high  using  the  consistency  argument.  It  appears  that  the  CCRC  pro¬ 
cedure  reduces  the  RMS .  However,  it  also  seems  to  add  to  the  bias 
error.  It  must  be  remembered  that  the  error  characteristics  (level  by 
level  independent  randomness)  used  are  the  severest  test  of  CCRC.  For 
a  more  trend  characteristic  error  as  opposed  to  random  error  the  CCRC 
works  better  as  evidenced  by  the  C02  values  in  Figure  34. 

In  conclusion,  this  analysis  has  shown  the  RMS  and  bias  errors  of 
Pw£  and  a  to  be  of  the  order  of  magnitude  given  in  Table  7.  These 
values  represent  acceptable  levels  when  matched  with  the  approach  of 
using  weighting  function  curves  defined  by  the  shape  parameters  Pwf  and 
a  as  an  indication  of  climatological  cloudiness .  As  pointed  out  above , 
a  great  advantage  is  gained  due  to  the  fact  that  maximum  curve  shape 
parameter  errors  occur  as  pw^  •+  1000  mb  and  as  a  -*■  0  and  for  these 
values  of  pw^  and  a  the  broadband  weighting  curve  shape  remains  nearly 
as  it  is  for  the  clear  sky  case. 

D.  Single  pw^  level  effect 

The  nature  of  this  error  effect,  due  to  errors  in  the  assumption 
that  the  satellite  sensor  is  viewing  a  scene  with  a  single  pw£  level, 
requires  a  different  analysis  approach.  Figure  38  illustrates  the 
effect  of  two  radiating  surfaces  on  the  spectral  weighting  function  for 


Spectral  weighting  curve  for  a  case  of  two  separate  radiating  surfaces 
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wavenumber  727.5  cm  ^  using  a  tropical  atmosphere.  The  values  of  Pwq 
and  Pw£2  are  210  and  510  mb  respectively.  The  a  values  are  0.35  and 
0.40  respectively.  It  is  clear  that  to  describe  the  shape  of  this 
weighting  function  profile  takes  4  curve  shape  parameters.  What  sort 
of  errors  are  obtained  when  the  set  of  spectral  radiance  values  that 
correspond  to  these  four  shape  parameters  are  used  in  the  2  parameter 
retrieval  program?  This  problem  of  two  radiating  surfaces  versus  one 
radiating  surface  is  addressed.  To  slightly  simplify  further,  a  values 
are  restricted  to  1.0  for  a  single  specified  known  scene  which  when 
paired  translates  into  a  specified  known  scene  that  is  described  by  a's 
of  0.5  and  0.5. 

As  pointed  out  in  Figure  1  the  ultimate  objective  is  to  composite 
broadband  weighting  curves  for  purposes  of  a  climate  index.  Specific 
scene  retrievals  are  only  a  secondary  consideration.  Therefore,  the 
particular  source  of  error  described  in  this  section  is  related  to  the 
final  composite  weighting  curve  product.  To  judge  error  effects,  the 
true  versus  retrieved  composite  curves  are  compared.  The  composite  can 
be  thought  of  as  representing  either  time  or  space  averages.  The  ex¬ 
periment  most  closely  resembles  a  space  average. 

To  do  this  comparison  the  following  experimental  procedure  is 
followed.  Assign  an  effective  radiating  surface  to  a  simulated  spot 
scene.  Do  this  for  many  spots  ('v  200  spots)  using  different  pw^  values 
with  a  *  1.  There  is  a  set  of  spectral  radiance  values  for  each  spot. 

A  composite  of  curves  200)  described  from  the  spot  scenes  spectral 
radiances  gives  the  true  weighting  curve.  Assume  the  satellite  sensor 
views  two  scenes  at  once.  In  other  words,  the  satellite  is  assumed  to 
view  two  effective  radiating  surfaces  instead  of  one.  Thus,  each  of 
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Che  true  spot  scenes  are  paired  off  (y  100  pairs) .  There  is  a  set  of 
spectral  radiances  for  each  pair  of  spot  scenes  which  represent  an  aver¬ 
age  of  the  two  separate  radiance  sets.  From  this  set  of  spectral  radi¬ 
ances,  pw£  and  a  values  are  retrieved  and  all  the  resulting  weighting 
curves  (n>  100)  are  composited  to  give  the  retrieved  composite  weighting 
function  curve.  This  is  done  for  three  separate  groups  of  200  spot 
values . 

Tables  16  and  17  give  the  effective  cloud  top  distributions  in  de¬ 
tail  for  the  three  groupings .  The  flat  distribution  (FLATS)  is  modeled 
after  20  day  average  GATE  B-array  cloud  data  (Cox  and  Griffith,  1978). 
The  other  two  distributions  (BIMOS  and  1PK35)  are  specified  to  resemble 
alternate  scenarios.  BIMOS  represents  a  specified  bimodal  distribution, 
while  1PK35  represents  a  specified  one  peak  distribution  at  the  200  mb 
layer.  There  are  5  levels  within  each  effective  cloud  top  category 
(i.e.  400,  420,  440,  460  and  480  mb).  The  spots  are  paired  off  as 
realistically  as  possible  by  category.  Within  each  pair  the  difference 
between  tops  vary  according  to  the  values  at  the  bottom  of  Table  8. 

For  example,  if  200  and  300  nb  category  tops  are  paired  off  together 
their  tops  may  differ  between  20  and  180  mb;  a  280  mb  top  may  be  paired 
with  a  300  mb  top  or  a  200  mb  top  may  be  paired  with  a  380  mb  top. 

Once  the  top  categories  are  chosen,  the  exact  top  differences  are 
determined  randomly.  The  exact  specified  distribution  of  top  differ¬ 
ence  values  is  largely  dependent  on  the  total  distribution  of  effective 
cloud  tops.  Notice  that  two  120  pair  models  of  the  flat  distribution 
are  specified  in  order  to  judge  result  sensitivity  to  differences  in 
pair  separation  distributions. 
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Table  16.  Model  effective  cloud  top  pair  separation  distribution. 
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Because  one  would  not  expect  each  and  every  set  of  satellite 
radiance  values  to  be  Influenced  by  a  scene  with  two  different  radiat¬ 
ing  surfaces,  this  test  set  up  may  be  extreme.  On  the  other  hand, 
occassional  scenes  with  more  than  two  radiating  surfaces  may  be  viewed. 
Overall  one  may  consider  this  test  of  the  single  representative  p  ^ 
value  error  to  be  not  quite  a  worse  case  situation. 

Tables  16  and  17  give  the  effective  cloud  top  retrieved  distribu¬ 
tions.  Table  18  shows  that  kind  of  single  pw^  and  a  values  the  retriev¬ 
al  program  calculates  from  sets  of  spectral  radiances  containing  the 

effects  of  two  p  ,  and  ot  values.  In  most  cases  the  results  seem  to 
*wf 

provide  a  fortuitous  average  type  value.  Figure  39  gives  a  more  in¬ 
formative  summary  of  this  process.  Plot  A  shows  the  three  effective 
cloud  top  distributions.  Plot  B  compares  the  specified  top  distribu¬ 
tion  (240  spots)  against  the  retrieved  top  distribution  (120  pairs). 
Notice  that  the  retrieved  top  distribution  fails  to  match  the  specified 
one  in  the  lower  troposphere.  This  fault  is  minimized  by  considering 
that  effective  clouds  at  that  level  have  little  influence  on  the  over¬ 
all  broadband  weighting  function.  Notice  that  both  pair  models  show 
the  same  general  results.  Plots  C  and  D  also  show  favorable  comparison 
between  the  distribution  of  specified  effective  cloud  tops  and  the  re¬ 
trieved  ones . 

Figures  40-42  show  the  broadband  weighting  curve  composites  for 
the  three  effective  cloud  top  distributions.  Specified  (i.e.  240  spot) 
curves  are  compared  to  retrieved  (i.e.  120  pair)  curves.  In  Figure  40 
the  120  pair  model  weighting  cuve  is  not  plotted  in  order  to  make 
visual  comparison  between  the  two  plotted  curves  easier.  All  three 
figures  show  only  minor  differences  between  the  specified  and  retrieved 
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Table  18.  Satellite  Spot  to  Pair  Value  Conversion  Example. 
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Figure  39.  Plots  showing  comparison  between  specified  effective  cloud 
top  distributions  and  the  respective  retrieved  distribu¬ 
tions  . 
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Broadband  weighting  curve8  for  specified  and  retrieved  flat  effective  cloud  top  distribution 
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Figure  41.  Broadband  weighting  curves  for  specified  and  retrieved  two  peaked  effective  cloud  top  distribu¬ 
tion. 
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Broadband  weighting  curves  for  specified  and  retrieved  one  peak  effective  cloud  top  distribu- 


curves.  Notice  that  errors  associated  with  u-per  tropospheric  radiat¬ 
ing  surface  retrievals  are  the  greatest.  Figure  41  allows  direct  com¬ 
parison  between  high  and  low  atmospheric  radiating  surfaces  and  their 
effect  on  the  weighting  curve.  Consequently  Figure  42,  which  is  the 
upper  tropospheric  single  peak  distribution  plot,  is  taken  to  represent 
the  severest  test  case.  If  retrieval  errors  are  acceptable  for  this 
case,  then  one  would  expect  them  to  be  acceptable  for  nearly  all  other 
effective  cloud  top  distributions. 

Table  8  gives  the  results  of  a  statistical  analysis  of  the  dif¬ 
ferences  between  the  plotted  curves  in  Figures  40-42.  The  statistical 
test  is  fully  explained  in  Chapter  VI.  The  values  in  Table  8  indicate 
that  for  all  cases  including  the  so-called  worse  case  there  is  a  high 
probability  that  the  specified  versus  retrieved  curves  are  not  dif¬ 
ferent. 

From  this  experiment  it  appears  that  the  assumption  of  one  repre¬ 
sentative  pwf  level  in  the  radiometer  field  of  view  leads  to  only  small 
errors  when  considered  for  the  case  of  compositing  weighting  function 
curves.  Of  course  for  the  single  retrieval  case,  calculating  a  single 
set  of  pw£  and  a  values  can  lead  to  a  retrieved  weighting  curve  that  is 
much  different  from  the  true  one.  There  may  be  a  question  of  the 
representativeness  of  the  effective  cloud  top  distributions  specified. 
Figure  16  gives  distribution  plots  of  window  radiances  which  for  the 
most  part  relate  directly  to  the  cloud  top  distribution.  Each  plot 
represents  over  900  measurements  from  a  200  x  200  km  section  for  a 
specific  time  from  a  geostationary  satellite.  Notice  the  large  number 
of  1  peak  distributions.  Although  errors  in  this  window  radiance  data 
analysis  tend  to  broaden  the  peaks,  they  are  still  quite  steep  peaks. 
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It  appears  from  an  examination  of  Figure  16  that  the  experimental  dis¬ 
tributions  specified  above  represent  nearly  a  worse  case  consideration. 

In  other  words,  steep  peaked  distributions  imply  that  the  assumption  of 
a  single  radiating  surface  is  for  many  cases  nearly  true.  The  tests 
above  address  the  problem  of  the  assumption  always  being  false. 

E.  Minimizing  the  errors 

Throughout  this  Appendix  different  techniques  have  been  used  to 
minimize  retrieval  weighting  function  errors  (sometimes  in  the  form  of 
pw£  and  a)  due  to  the  shortcomings  of  the  basic  approach  assumptions. 

These  methods  include  using  a  clear  column  radiance  correction  (CCRC) 
and  stressing  the  use  of  composite  weighting  curves  instead  of  single 
scene  retrievals.  The  treatment  of  boundary  conditions  explained  in 
the  Chapter  V  introduction  may  also  be  viewed  as  a  device  to  limit  re¬ 
trieval  errors. 

It  is  also  important  to  view  the  general  approach  of  using  weight¬ 
ing  function  curves  as  a  major  way  of  minimizing  errors.  Used  as  a 
representation  of  the  atmosphere's  vertical  infrared  radiative  emitting 
structure  (VIRES),  the  composite  broadband  weighting  curves  are  a  meas¬ 
ure  of  climatological  cloudiness  that  is  most  complimentary  to  the 
satellite  radiative  measurements.  As  a  result,  the  conditions  under 
which  the  retrieval  of  weighting  function  curve  shape  parameters  is 
weakest  is  exactly  the  condition  under  which  the  error  effect  is  least. 

An  advantage  important  to  the  approach  of  averaging  over  time  is  the 
ability  to  use  this  retrieval  technique  day  or  night  with  equal  reli¬ 
ability. 
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Finally,  to  minimize  errors  further  one  uses  the  best  Information 
possible.  This  Includes  using  retrieved  clear  sky  radiance  values.  A 
disadvantage  of  this  technique  is  that  to  use  retrieved  temperatures 
and  humidity  profiles  would  mean  solving  for  the  C  and  D  values  again. 

Of  course  sets  of  C  and  D  values  for  different  scenarios  could  be  pre¬ 
calculated.  Then  the  set  corresponding  closest  to  the  measured  pro¬ 
files  could  be  used.  However,  for  most  applications  simple  season 
values  of  C  and  D  for  specific  areas  should  be  sufficient  when  used  with 
CCRC. 

In  suranary,  this  Appendix  examined  the  errors  associated  with  re¬ 
trieval  of  the  weighting  curve  shape  parameters  or  with  the  weighting 
curve  composites.  As  indicated  by  Tables  7  and  8  the  empirical  RTE 
technique  gives  results  that  in  the  mean  represent  small  errors.  Of 
course  only  a  sample  of  reasonable  simulated  error  source  values  are 
used  in  the  analysis.  Nevertheless,  the  results  indicate  that  the 
empirical  RTE  retrieval  technique  gives  usefully  accurate  results. 


